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Abstract 

Modified Newtonian dynamics (MOND) is a possible way to explain the flat galaxy rotation 
curves without invoking the existence of dark matter. It is however quite difficult to predict such a 
phenomenology in a consistent field theory, free of instabilities and admitting a well-posed Cauchy 
problem. We examine critically various proposals of the literature, and underline their successes and 
failures both from the experimental and the field-theoretical viewpoints. We exhibit new difficulties 
in both cases, and point out the hidden fine tuning of some models. On the other hand, we show 
that several published no-go theorems are based on hypotheses which may be unnecessary, so that 
the space of possible models is a priori larger. We examine a new route to reproduce the MOND 
physics, in which the field equations are particularly simple outside matter. However, the analysis 
of the field equations within matter (a crucial point which is often forgotten in the literature) 
exhibits a deadly problem, namely that they do not remain always hyperbolic. Incidentally, we 
prove that the same theoretical framework provides a stable and well-posed model able to reproduce 
the Pioneer anomaly without spoiling any of the precision tests of general relativity. Our conclusion 
is that all MOND-like models proposed in the literature, including the new ones examined in this 
paper, present serious difficulties: Not only they are unnaturally fine tuned, but they also fail to 
reproduce some experimental facts or are unstable or inconsistent as field theories. However, some 
frameworks, notably the tensor-vector-scalar (TeVeS) one of Bekenstein and Sanders, seem more 
promising than others, and our discussion underlines in which directions one should try to improve 
them. 

PACS numbers: 04.50. +h, 95.30.Sf, 95.35.+d 
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I. INTRODUCTION 



Although general relativity (GR) passes all precision tests with flying colors [H, 0], there 
remain some puzzling experimental issues, notably the fact that the Universe seems to be 
fllled with about 72% of dark energy and 24% of dark matter, only 4% of its energy content 
being made of ordinary baryonic matter [3|,0]- Dark energy, a fluid whose negative pressure is 
close to the opposite of its energy density, would be responsible for the present accelerated 
expansion of the Universe [1, 0, 0, @|- Dark matter, a pressureless and noninteracting 
component of matter detected only by its gravitational influence, is suggested by several 
experimental data, notably the flat rotation curves of clusters and galaxies Another 
experimental issue is the anomalous extra acceleration 5a ~ 8.5 x 10~^°m.s~^ towards the 



Sun that the two Pioneer spacecrafts exhibited between 30 and 70 AU [lOj (see also 

Actually, none of these issues contradicts GR nor Newtonian gravity directly. Indeed, 
dark energy may be understood as the existence of a tiny cosmological constant A 3 x 
10~^^^ c" lihG), or as a scalar fleld (called quintessence) slowly rolling down a potential 
T^ . 13 . Many candidates of dark matter particles have also been predicted by different 



theoretical models (notably the class of neutralinos, either light or more massive, occurring 



in super symmetric theories; see e.g. [I5j), and numerical simulations of structure formation 



have obtained great successes while incorporating such a dark matter (see, e.g., [l6|). [It 
remains however to explain why A is so small but nonzero, and why the dark energy and 
dark matter densities happen to be today of the same order of magnitude.] The Pioneer 
anomaly seems more problematic, but the two spacecrafts are identical and were not built 
to test gravity, therefore one must keep cautious in interpreting their data. A dedicated 
mission would be necessary to conflrm the existence of such an anomalous acceleration. 

Nevertheless, these various issues, considered simultaneously, give us a hint that Newton's 
law might need to be modifled at large distances, instead of invoking the existence of several 
dark fluids. To avoid the dark matter hypothesis, Milgrom proposed in 1983 such 
a phenomenological modiflcation, which superbly accounts for galaxy rotation curves [l8| 
(although galaxy clusters anyway require some amount of dark matter), and automatically 
recovers the Tully-Fisher law [19| f ^ oc M (where M denotes the baryonic mass of a galaxy, 
and foo the asymptotic circular velocity of visible matter in its outer region). The norm 
a of a particle's acceleration is assumed to be given by its Newtonian value a^v when it is 
greater than a universal constant Oq ~ 1.2 x 10~^°m.s~^, but to read a = ya/vOo in the 
small-acceleration regime a < qq. In particular, the gravitational acceleration should now 
read a = \/GMa^/r at large distances, instead of the usual GM/r"^ law. 

Various attempts have been made to derive such a modifled Newtonian dynamics 
(MOND) from a consistent relativistic fleld theory. The main aim of the present paper 
is to examine them critically, by underlining their generic difficulties and how some models 
managed to solve them. Three classes of difficulties may actually be distinguished: (i) theo- 
retical ones, namely whether the proposed model derives from an action principle, is stable 
and admits a well-posed Cauchy problem; (ii) experimental ones, in particular whether 
solar-system and binary-pulsar tests are passed, and whether the predicted light deflection 
by galaxies or haloes is consistent with weak-lensing observations; and (iii) "esthetical" ones, 
i.e., whether the proposed model is natural enough to be considered as predictive, or so flne- 
tuned that it almost becomes a flt of experimental data. These different classes of difficulties 
are anyway related, because fine-tuning is often necessary to avoid an experimental problem, 
and because some tunings sometimes hide serious theoretical inconsistencies. 
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We will not address below the problem of dark energy, although it may also be tackled 



with a modified-gravity viewpoint. For instance the DGP brane model [20|, |21[ predicts an 
accelerated expansion of the Universe without invoking a cosmological constant 22|, |23|], and 
the k-essence models 2J, |25|, |26| also reproduce many cosmological features in a somewhat 
more natural way than quintessence theories. Some of our theoretical discussions below, 
notably about stability and causality, are nevertheless directly relevant to such models. In- 
deed, k-essence theories are characterized by a non-linear function of a scalar field's kinetic 
term, precisely in the same way relativistic aquadraric Lagrangians (RAQUAL) were de- 
vised to reproduce the MOND phenomenology in a consistent field theory [27|. Our aim is 
not to discuss the Pioneer anomaly in depth either. Although the numerical value of these 
spacecrafts' extra acceleration 5a is of the same order of magnitude as the MOND parameter 
oq, there indeed exist important differences between their behavior and the MOND dynam- 
ics. However, we will come back to the Pioneer anomaly at the very end of the present 
paper. While analyzing a new class of models a priori devised to reproduce the MOND 
phenomenology, we will show that it can account for the Pioneer anomaly without spoiling 
any of the precision tests of general relativity, while being stable and admitting a well-posed 
Cauchy problem. 

Our paper is organized as follows. In Sec. [Tll we examine the various ways one may try to 
modify the laws of gravity in a consistent relativistic field theory, and give a critical review 
of several models proposed in the literature. We notably underline that an action should not 
depend on the mass of a galaxy, otherwise one is defining a different theory for each galaxy. 
We also recall why higher-order gravity is generically unstable, contrary to scalar-tensor 
theories. When the dynamics of the scalar field is defined by an aquadratic kinetic term 
(RAQUAL or k-essence theories), we discuss the different consistency conditions it must 
satisfy. We point out that such conditions suffice for local causality to be satisfied, although 
some modes may propagate faster than light or gravitational waves. We finally recall how 
Ref. [23] reproduced the MOND gravitational force thanks to a RAQUAL model, but we 
exhibit a serious fine-tuning problem which does not seem to have been discussed before — 
or at least was not pointed out so clearly. 

In Sec. we recall that many models fail to reproduce the observed light deflection by 
"dark matter" haloes. We exhibit a counter-example to an erroneous claim in the literature, 
although this counter-example does not reproduce the correct MOND phenomenology. We 
also recall how a "disformal" (i.e., non-conformal) coupling of the scalar field to matter 



allowed Refs. [28|, l29j to predict the right light deflection, but that this class of models has 
been discarded too quickly, because of the existence of superluminal gravitons. Actually, the 
consistency and causality of such models is clear when analyzed in the Einstein frame. On 
the other hand, we show that the consistency of the field equations within matter implies 
new conditions which must be satisfied by the functions defining the theory. As far as we 
are aware, these crucial conditions were not derived previously in the literature. We finally 



discuss how the best present model, called TeVeS [30|, iMl, l32, l33|, solved the problem of 
light deflection by considering a stratified theory, involving a unit vector field in addition to 
a metric tensor and one or several scalar fields. 

In Sec. llVt we discuss the various difficulties which anyway remain in this TeVeS model, 
including some which have already been discussed in the literature and that we merely 
summarize. Reference 34| proved that the TeVeS Hamiltonian is not bounded by below, 
and therefore that the model is unstable. We also point out other instabilities in related 
models. Many papers underlined that stratified theories define a preferred frame, and that 
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they are a priori inconsistent with solar-system tests of local Lorentz invariance of gravity. 
However, the vector field is assumed to be dynamical in TeVeS, and it has been argued in 
'm that such tests are now passed. Although we do not perform a full analysis ourselves, 



our discussion suggests that preferred-frame effects are actually expected in TeVeS, and 
that they may be avoided only at the expense of an unnatural fine-tuning. The reason why 
disformal models were discarded in favor of the stratified TeVeS theory was mainly to avoid 
superluminal gravitons. However, the opposite phenomenon occurs in TeVeS: photons (and 
high-energy matter particles) propagate faster than gravitons. References jssl . 36| actually 
proved that the observation of high-energy cosmic rays imposes tight constraints on such 
a behavior, because these rays should have lost their energy by Cherenkov radiation of 
gravitational waves. After recalling this argument, we conclude that the sign of one of the 
TeVeS parameters should be flipped, implying that gravitons now propagate faster than 
light, like in the RAQUAL models. We finally discuss briefly binary-pulsar constraints, still 
without performing a full analysis, but underlining that a large amount of energy should 
be emitted as dipolar waves of the scalar field. The matter-scalar coupling constant should 
thus be small enough in TeVeS to pass binary-pulsar tests, and this implies an unnatural 
fine-tuning of its Lagrangian. 

Section |V] is devoted to new models whose point of view differs significantly from those 
of the literature. They are extremely simple in vacuum, reducing respectively to general 
relativity and to Brans-Dicke theory. The MOND phenomenology is obtained via a non- 
minimal coupling of matter to the gravitational field(s). However, the first model exhibits a 
subtle instability, similar to the one occurring in higher-order gravity (although no tachyon 
nor ghost degree of freedom may be identified in vacuum). The second model avoids this 
instability due to higher derivatives, and reproduces the MOND phenomenology (including 
light deflection) quite simply, as compared to the literature. However, our analysis of the 
field equations shows that they do not remain hyperbolic within the dilute gas in outer 
regions of a galaxy. Therefore, although promising, this model is inconsistent. But the same 
framework provides a simple model to reproduce consistently the Pioneer anomaly, without 
spoiling the precision tests of general relativity. 

Finally, we present our conclusions in Sec. IVI[ We notably mention other experimental 
constraints that any complete MOND-like theory of gravity should also satisfy, besides those 
which are discussed in the present paper. 



II. LOOKING FOR MOND-LIKE FIELD THEORIES 



A. General relativity 

Einstein's general relativity is based on two independent hypotheses, which are most 
conveniently described by decomposing its action as S' = ^gravity + •S'matter- First, it assumes 
that all matter fields, including gauge bosons, are minimally coupled to a single metric tensor, 
that we will denote ^^j/ throughout the present paper. This metric defines the lengths and 
times measured by laboratory rods and clocks (made of matter), and is thereby called the 
"physical metric" (the name "Jordan metric" is also often used in the literature) . The matter 
action may thus be written as SmatteAi^', dfiu], where ip denotes globally all matter fields, 
and where the angular brackets indicate a functional dependence. This so-called "metric 
coupling" implies the weak equivalence principle, i.e., the fact that laboratory-size objects 
fall with the same acceleration in an external gravitational field, which is experimentally 
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verified to a few parts in 10^^ 37], |38|. For instance, tlie action of a point-particle reads 
S'pp = — / mc\J —g^y{x)v^^v'' dt, and depends both on the spacetime position of the 
particle and its velocity = dx^/dt. The second building block of GR is the Einstein- 
Hilbert action 

C d'^x 



5'gravity — , „ ^ / \/~9* R* i (2-1) 

ibvrG J c 

which defines the dynamics of a spin-2 field g*^^, called the "Einstein metric" . We use the sign 
conventions of [s^, notably the mostly-plus signature, and always indicate with a tilde or 
a star (either upper or lower) which metric is used to construct the corresponding quantity. 
For instance, g^, = det{g^^) is the determinant of the Einstein metric and R* its scalar 
curvature. Similarly, we will denote as and V* the covariant derivatives corresponding 
respectively to the Jordan and Einstein metrics, and as g^'^ and gf^'^ the inverses of these 
two metrics. This rather heavy notation will allow us to be always sure of what we are 
talking about in the following. Einstein's second hypothesis is that both metrics coincide: 

9fMu = gf_iu- 



Milgrom analyzed in [40|, |41j the consequences of modifying the matter action, that he 
called "modified inertia" . He focused on a point particle in an external gravitational field, 
and assumed that its action could depend also on the acceleration & = d\/dt and its higher 
time-derivatives. However, he proved that to obtain both the Newtonian and the MOND 
limits and satisfy Galileo invariance, the action must depend on all time-derivatives d'^v/dt'^ 
to any order, i.e., that the action is necessarily nonlocal. This does not necessarily violate 
causality (see the counter-example in ^421]), and is actually a good feature for the stability 
of the theory (see Sec. Ill CI below), but the actual computation of the predictions is quite 
involved. We refer the reader to the detailed paper [401] for more information about this 
interesting viewpoint, but we focus below on metric field theories, i.e., such that matter is 
minimally coupled to g^u as in GR, with only first derivatives of the matter fields entering 
the action. 

Another way to modify GR is to assume that this physical metric does not propagate as a 
pure spin-2 field, i.e., that its dynamics is no longer described by the Einstein- Hilbert action 
f l2.ip . and that g^i, is actually a combination of various fields. Since the weak equivalence 
principle is very well tested experimentally, and implied by a metric coupling, most MOND- 
like models in the literature focused on such a "modified gravity" viewpoint, and we will 
examine them in the present paper. The fact that they involve extra fields, besides the usual 
graviton (fluctuation of g*^^,) and the various matter fields entering the Standard Model of 
particle physics, makes the distinction with dark matter models rather subtle. The crucial 
difference is that, in the dark- matter paradigm, the amount of dark matter is imposed by 
initial conditions, and its clustering generates gravitational wells in which baryonic matter 
fall to form galaxies and large-scale structures. On the other hand, in the modified-gravity 
viewpoint, baryonic matter generates itself an effective dark-matter halo Mdark oc ■>/ Mbaryon- 
Such a halo may just be an artifact of the way we interpret the gravitational field of baryonic 
matter alone at large distances. But it may also be a real dark-matter halo, made of the extra 
gravitational fields, and generating itself a Newtonian potential. In such a case, the difference 
with standard dark-matter models would be that its mass Mdark is imposed as above by the 
baryonic one, and modified gravity could thus be considered as a constrained class of dark- 
matter models. Most of the theories that we will mention in the following predict that the 



energy density of the extra gravitational degrees of freedom is negligible with respect to 

pextra fields I ^ iT^ma 



baryonic matter: [j'extra fields | ^ jjimatterj^ Therefore, the effective dark-matter haloes will be 



5 



most of the time artifacts of our interpretation of observations. However, we will also consider 
in Sec. [mB] a model of the "constrained dark matter" type (|T;^^t>-a fields | ^ |j.matter|^^ ^ 

counter-example to some claims in the literature about light deflection.^ 

As suggested by Milgrom, one may also consider models modifying both inertia and 
gravity. Paradoxically, as we will see in Sec. IV Al below, there also exists a non-trivial (non- 
GR) possibility in which neither of them is modified in the above sense. Matter will be 
universally coupled to a second-rank symmetric tensor g^^, i.e., still described by an action 
'S'matter ['?/'; ^^j!/], and the dynamics of gravity will be described by the pure Einstein- Hilbert 
action (12.11) in vacuum, but the two metrics and (7*^ will nevertheless differ within matter. 



B. Various ideas in the literature 



Although interesting from a phenomenological point of view, some models proposed in 
the literature write field equations which do not derive from an action, and cannot be 
obtained within a consistent field theory. For instance, actually not to reproduce a MOND- 
like behavior but anyway as a model of dark matter, Ref. proposes to couple differently 
dark matter and baryonic matter to grav ity (an idea explored in the cosmological context by 
various authors, notably (iT 



47, 48.1491. [sol. Isil. 52, S^]). The authors underline themselves 
that it cannot be considered as a fundamental theory (and they actually need some negative 
energy to obtain a repulsive force), but it is also instructive to stress that it cannot derive 
from an action. Indeed, the scalar field entering this model is assumed not to be generated 
by baryons, which implies that the scalar-baryon vertex vanishes. On the other hand, the 
assumed equation of motion for this baryonic matter does depend on the scalar field, thereby 
implying that the scalar-baryon vertex does not vanish. In conclusion, although the behavior 
of such a model might be mimicked as an effective theory of a more fundamental one, it 
cannot be described itself as a field theory, and notably does not satisfy conservation laws. 
In the present paper, our discussion will be restricted to consistent field theories. 

Some models in the literature do write actions, but they dep end on the (baryonic or 



luminous) mass M of the galaxy [5^, |55j. For instance Ref. [55|] proposes a gravitational 
action in which the kinetic term is some power of the Ricci scalar Since this implies, 

in general, an asymptotic velocity of the form v'^ oc nc^ + 0{GM/r), Ref. 55|] then proposes 



to replace n by a/ GMao / in order to recover the Tully-Fisher law. However, without an 
additional prescription (which must be non-local in nature) specifying which mass M enters 
the action, this theory is ill defined. We will come back to such f{R) theories of gravity in 
Sec. Ill CI and in Sec. IIIDI It should be also underlined that a gravitational force oc 1/r is 
not difficult to obtain by choosing an appropriate scalar field potential, for instance, as we 
will illustrate in Sec. HID] below. However, the second crucial feature of MOND, which arises 
from the Tully-Fisher law, is also the factor y/GMa^ oc \fM multiplying this 1/r. Some 
consistent field theories actually obtain a flattening of rotation curves but do not predict 



the ri ght amplitude of the asymptotic velocity [56|, l57 



In |58l . l59| , Einstein and Straus studied an extension of general relativity in which the met- 



ric tensor is nonsymmetric, g^^ 7^ g^^. In modern language, this corresponds to considering 



^ Let us also mention that the recent reinterpretation of MOND proposed in Refs. 43], |4J| is also of the 
constrained dark matter type: The MOND force oc 1/r is caused by a background of gravitational dipoles 
contributing directly to the Newtonian potential. 
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an antisymmetric tensor field as a partner to the usual graviton g{fj,u)- Moffat analyzed 
the phenomenology of this class of models in many papers (starting with [60|l), and it was 
shown to define a consistent and stable field theory provided B^^^j is massive |6l| . More re- 
cently, Moffat showed that it may reproduce the fiat rotation curves of clusters and galaxies 



[62, 1631, [64, l65|, l66|, l67j, probably even better that the original MOND proposal. However, 
several assumptions are needed to obtain such a prediction, and it cannot yet be considered 
as a predictive field theory, in our opinion. Indeed, the author points out himself that a 
constant entering his action must take three different values at the solar system, galaxy, and 
cluster scales. It must therefore be considered as a running (renormalized) coupling constant 
instead of a pure number imposed in the action, and a complete theory should be able to 
predict the three needed values. A related difficulty is the fact that the force generated 
by the antisymmetric tensor field B^^^^ is repulsive, whereas galaxy rotation curves need a 
gravitational attraction greater that the Newtonian one. Here again, the author invokes 
renormalization of the constants to obtain the right behavior, but does not provide a full 



derivation. But the most instructive difficulty of Ref. [63j is that the extra force is predicted 
to be cx kM'^/r instead of the needed cx y/M/r behavior. The author first assumed in [g^] 
that an unknown mechanism could tune the proportionality constant to be A; oc M~^/^, so 
that the right MOND phenomenology would be recovered. More recently [g^, he managed 
to derive this relation by assuming a particular form of a potential. However, this potential 
does depend explicitly on the mass M of the galaxy, and the predictivity of the model is 
thus lost. In conclusion, although the framework of nonsymmetric gravity seems promising, 
it is not yet formulated as a closed consistent field theory, and it should only be considered 
at present as a phenomenological fit of observational data. 

In the following, we will focus on field theories whose actions do not depend on the mass 
or scale of the considered objects. They may involve fine-tuned numbers, but they will be 
fixed once for all. This is necessary to define a predictive model, but we will underline that 
it does not suffice to prove its consistency. In particular, the stability needs to be analyzed 
carefully. 



C. Higher-order gravity 

One-loop divergences of quantized GR [68|] are well known to generate terms proportional 
to^ R^, R^^ and R^^yp^j- It is thus natural to consider extensions of GR already involving 
such higher powers of the curvature tensor in the classical action. At quadratic order, it is 
always possible to write them as 

^gravity = I^V^[R + a Cl^p^ + /3 + 7 GB] , (2.2) 

where C^ypc, denotes the Weyl (fully traceless) tensor, GB = R^^^p^ — 4i?^^, + R^ is the 
Gauss-Bonnet topological invariant, and 7 are constants having the dimension of a 
length squared (i.e., (h/mc)'^ if m denotes the corresponding mass scale). The topological 
invariant GB does not contribute to the local field equations and may thus be discarded. 
In his famous thesis j69|, Stelle proved that such an action gives a renormalizable quantum 



^ To simplify, we do not write the tildes which should decorate all quantities in this subsection, indicating 
that 3^1/ = 3^1/ will be later assumed to be minimally coupled to matter. 
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theory, to all orders, provided both a and /3 are nonzero. However, he also underlined 
that this cannot be the ultimate answer to quantum gravity because such a theory contains 
a ghost, i.e., a degree of freedom whose kinetic energy is negative. To understand this 
intuitively, it suffices to decompose the schematic propagator in irreducible fractions: 

111 , ^ 

(2.3) 



where the first term, 1/p^, corresponds to the propagator of the usual massless graviton, 
whereas the second term corresponds to a massive degree of freedom such that m? = 1/a. 
The negative sign of this second term indicates that it carries negative energy. Note that 
one cannot change its sign by playing with that of a: If a < 0, this extra degree of freedom 
is even a tachyon (negative mass squared), but it anyway remains a ghost (negative kinetic 
energy). The problem with such a ghost degree of freedom is that the theory is violently 
unstable. Indeed, the vacuum can disintegrate in an arbitrary amount of positive-energy 
usual gravitons whose energy is balanced by negative-energy ghosts. Even at the classical 
level, although such an instability is difficult to exhibit explicitly on a toy model, one expects 
anyway the perturbations of a given background to generically diverge, by creating growing 
gravitational waves containing both positive-energy (massless) modes and negative-energy 
(massive) ones. 

The actual calculation, taking into account all contracted indices in the propagators, 
confirms the above schematic reasoning for the cxCj^ypa term in action (12.21) . and shows that 
the extra massive mode has also a spin 2, like the usual massless graviton. On the other 
hand, the (3 term generates in fact a positive-energy massive scalar degree of freedom. We 
will recall a well-known and simple derivation in Sec. Ill Dl below. The intuitive explanation 
is that the scalar mode corresponds in GR to the negative gravitational binding energy 
(which is not an actual degree of freedom because it is constrained by the field equations), 
so that the negative sign entering the right-hand side of Eq. (12.31) in fact multiplies an already 
negative term. The resulting extra degree of freedom thereby appears as a positive-energy 
scalar field 



More general actions /(-R, R^u, R^upa) have been considered several times in the literature. 



notably in [7l| and [72[. The conclusion is that they generically contain the massive spin-2 
ghost exhibited above, which ruins the stability of the theory. The only allowed models, 
within this class, are functions of the scalar curvature alone, /(-R), that we will study in 
Sec. [TID] below. 

The fact that such functions of the curvature generically involve a ghost tells us that some 
particular models may avoid it. One may for instance consider a Lagrangian of the form 
i?-|-7GB + a;(-R^j^p^)", n >2, whose quadratic term reduces to the Gauss-Bonnet topological 
invariant GB = R^^yp^j^'^R^^u + R^ ■ Around a fiat Minkowski background, the only quadratic 
kinetic term 0{gp,i, — rip^Y is thus the one coming from the standard Einstein-Hilbert term 



i?, so that no propagator can be defined for any ghost degree of freedom. However, Ref. [71 
underlined that fiat spacetime is generically not a solution of such higher-order gravity 
theories. Around a curved background, the second-order expansion of higher-order terms 
like C^L-pcr)" thereby generates a nonzero kinetic term for the spin-2 ghost (see Fig. [H] in 
Sec. IV Al below for a diagrammatic illustration). 

One may still try to devise higher-order gravity models such that their second-order 
expansion around any background never generates a negative-energy kinetic term. For 



instance, Refs. [73|, [T^, l75| consider Lagrangians of the form R + / (GB), and show that the 
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two degrees of freedom corresponding to the generic spin-2 ghost are not excited. A similar 
trick as the one presented in Sec. Ill Dl below for f{R) models indeed proves that R + /(GB) 
theories only involve a scalar degree of freedom in addition to the usual massless spin- 
2 graviton. However, the kinetic terms take nonstandard expressions in this framework, 
and at present there is no proof that the full Hamiltonian is bounded by below. One 
can even show that ghost modes do exist in particular backgrounds (see, e.g., Eq. (22) of 
Ref. [tsI), although this does not prove either that the Hamiltonian is necessarily unbounded 
by below. In our opinion, the clearest hint that R + /(GB) models are probably unstable 
comes from the schematic reasoning of Eq. (12. 3p above. In f{R) models, the scalar degree of 
freedom happened to have positive kinetic energy because it corresponded to a higher-order 
excitation of the negative-energy Newtonian potential (constrained by the field equations). 
In R + /(GB) models, the scalar degree of freedom a priori corresponds to a higher-order 
excitation of another mode present in GR, generically dynamical and therefore carrying 
positive energy. The crucial minus sign entering Eq. (12. 3p then shows that the new scalar 
mode should probably be a ghost. 

More generally, deadly instabilities can exist in higher-order field theories even if no ghost 
degree of freedom can be identified in a perturbative way. Indeed, as recalled notably in 
their Hamiltonian is generically unbounded by below because it is linear in at least 



one of the canonical momenta. As in [70j, let us illustrate this on the simplest example of 
a Lagrangian depending on a variable q and its first two time derivatives, say C{q,q,q). If 
q can be eliminated from C by partial integration, we are in the standard case of a theory 
depending only on first derivatives, and we know that some models do give a Hamiltonian 
which is bounded by below. This is the case when q appears only linearly in C, even 
multiplied by a function of q and q since J dt q'^cf^q = —nj (m + 1) J dt g"-ig»^+2_|_ boundary 
terms. Let us thus consider only the case of a "non-degenerate" Lagrangian,^ i.e., such that 
the definition p2 = dC/dq can be inverted to express q as a function of g, q and p2'- 

g = /(g,g,P2). (2.4) 



In such a case, Ostrogradski showed in 1850 [76[ that the initial data must be specified by 
two pairs of conjugate momenta, that he defined as 

_dC d (dC\ 

^^ = ^' ^' = Tq-jA-d-q)^ ^'-'^^ 
dC 

q2 = q, p2 = -^, (2.5b) 

and he proved that the following definition of the Hamiltonian 

n = piqi + P2q2 - C{q, q, q) (2.6) 

does generate time translations. Indeed, qi = dH/dpi and pi = —dTC/dqi reproduce the 
Euler-Lagrange equations of motion deriving from the original Lagrangian C However, this 
Hamiltonian must be expressed in terms of the momenta q^ and Pi defined in Eqs. (12. 5p 
above. Recalling that Eq. (12.40 allows us to write q = f{qi,q2,P2), one gets 

^ = Piq2+P2f{quq2,P2) - ^{quq2,f{qi,q2,P2)) ■ (2.7) 



^ Note that the R + /(GB) models studied in Refs. are degenerate, so that their Hamiltonian 

cannot be written in the Ostrogradski form (|2.6p . but this does not suffice to prove their stability. 
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The crucial problem of this expression is that it is linear in pi, and therefore unbounded by 
below [tO]. In other words, the theory is necessarily unstable, even if one does not identify 
any explicit contribution —pf defining a ghost degree of freedom. If one diagonalizes the 
kinetic terms in Eq. (12. 7p at quadratic order, say in terms of new momenta p^, then the 
standard positive-energy degree of freedom should appear as +pf . Since we know that H 
is unbounded by below, the other momentum must have a negative contribution, but it 
may be for instance of the form —p'2, and therefore absent at quadratic order. 

This result can be straightforwardly extended to (non-degenerate) models depending on 
even higher time-derivatives of q, say up to d^q/dt^. In that case, Ref. j?^ shows that 
the Hamiltonian is linear in n — 1 of the momenta pi, and thereby unbounded by below. 
On the other hand, Ostrogradski's construction of the Hamiltonian cannot be used if C 
depends on an infinite number of time derivatives, i.e., if it defines a nonlocal theory. In 
such a case, the theory may actually be stable although its expansion looks pathological [t^] ■ 
This is notably the case in the effective low-energy models defined by string theory, which 
do involve quadratic curvature terms like Eq. fl2.2l) above, but also any higher derivative 
of the curvature tensor (whose phenomenological effects occur at the same order as the 
quadratic terms). Nonlocal models of MOND have been studied in the literature [i^, and 
they can be proved to satisfy anyway causality, but they remain difficult to study from a 
phenomenological point of view. In the following, we will focus on local field theories. 



D. Scalar-tensor theories 



As mentioned in the previous section, gravity models whose Lagrangians are given by 
functions of the scalar curvature, f{R), do not exhibit any ghost degree of freedom, and 
avoid the generic "Ostrogradskian" instability of hi ghe r-derivative theories [70|. Let us 



recall here the simplest way to prove it (see 




and note 



that the following derivation assumes that the scalar curvature i? is a function of the metric 
Qfj^i, and its derivatives alone; the derivation and the result differ in the first-order Palatini 
formalism, where the scalar field does not acquire any kinetic term 85|, l86|). We start from 
an action 

d^x 



gravity 



gfiR), 



(2.8) 



where the global factor c^/lGvrG is temporarily set to 1, 
Lagrange parameter to rewrite this action as 



to simplify. We now introduce a 



S, 



gravity 



d X 



/(0)+ /(0) 



(2.9) 



The field equation for reads — 0j f"{4>) = 0, and implies (j) = R within each spacetime 

domain where /"(</>) 7^ 0. [If there exist hypersurfaces or more general domains where 
/"(0) = 0, then either / is constant and there is no gravitational degree of freedom in such 
domains, or f{R) oc i? — 2A and the theory reduces to GR plus a possible cosmo logical 
constant. In both cases, the scalar degree of freedom that we will exhibit below cannot 
propagate within such domains. A consistent field theory, without any discontinuity of its 
degrees of freedom, can thus be defined only within the domains, possibly infinite, where 
f"{R) 7^ 0.] The field equations for the metric g^i, may now be derived from action (12.91) . 
and one finds that they reduce to those deriving from action (12. 8p when is replaced by R. 
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This can also be seen by the a priori ilhcit use of the field equation cf) = R directly within 
action (12. 9p . which reduces trivially to (12.80 . In other words, the theory defined by action 
(12. 8p is equivalent to the scalar-tensor one 

^gravity = j d'x ^ - (9^0)' - [0/(0) - /(0)]} , (2.10) 

within each spacetime domain where f"{R) never vanishes. If one redefines the scalar field 
as $ = fifp), this action takes the form of the famous (Jordan-Fierz)-Brans-Dicke theory 
871 . l88l . |89| with no explicit kinetic term for i.e., with a vanishing ubd parameter, and 
with a potential defined by the last term within square brackets. It should be stressed 



that solar-system experiments [90| impose the bound cubd > 4000 when the scalar-field 
potential vanishes. Therefore, if matter is assumed to be coupled to g^^, the resulting 
scalar degree of freedom needs to be massive enough to have a negligible effect in the solar- 
system gravitational physics. The shape of the initial function f{R) is thus constrained by 
experiment. 

Although action (12.101) does not involve any explicit kinetic term for (or the Brans-Dicke 
scalar $), this scalar field does propagate. Indeed, the first contribution f'{(f))R contains 
terms of the form (pd'^g, i.e., cross terms —d(f)dg after partial integration. A redefinition of 
the fields then allows us to diagonalize their kinetic terms. This is achieved with the new 
variables 

g;^ = fm,., (2.11a) 

^ ^ ^ln/'(0), (2.11b) 

, , 0f (0) - f(0) 

and the full action, Eq. (I2.10p plus the matter part Smatteri^; ^^i/], now reads 

= 4^/ ^^{t " l^*'^^^^^'^^ - ^(^)} + ^mattej^;^,. = A\^)glX (2-12) 

where we have put back the global factor c" /IQ-kG multiplying the gravitational action, 
and where A{ip) = e'^^^. This is a particular case of scalar-tensor theory [9l|, where the 
precise matter-scalar coupling function A{(p) = e'fl^ comes from our initial hypothesis of 
a /(-R) theory. More general models may involve an arbitrary nonvanishing function A (99) 
relating the Einstein metric g*^^, (whose fiuctuations define the spin-2 degree of freedom) 
to the physical one g^^_y = A'^{ip)g*^^. Action (I2.12p clearly shows that the spin-0 degree of 
freedom ip does propagate and carries positive kinetic energy. Of course, this does not suffice 
to guarantee the stability of the theory: The potential V{(p), Eq. (I2.11cp . also needs to be 
bounded by below, and this imposes constraints on the initial function f{R)- 

Although such tensor-mono-scalar models are perfectly well-defined field theories, they 
do not reproduce the MOND phenomenology, at least not in the most natural cases. Indeed, 
when the potential V{ip) has a negligible infiuence, the scalar-field contribution to the grav- 
itational interaction is proportional to the Newtonian one, at lowest order, i.e., of the form 
(f oc GM/rc^ where M is the baryonic mass of the body generating the field. Therefore, 
it does not give the MOND potential 1/ GMao In r we are looking for. It has been shown 
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in 92, 93, 93] that Gaussian matter-scalar coupling functions A'^{ip) = e^°'^^ with Pq < 0, 
generate nonperturbative effects such that if is no longer strictly proportional to M, but 
still almost so for two different mass ranges. Therefore, the factor a/M we are looking for 
cannot be obtained that way either, besides the fact that the radial dependence of the scalar 
field is still oc 1/r instead of being logarithmic. When the potential V{(p) has a minimum, 
its second derivative defines the mass m of the scalar degree of freedom, and its contri- 
bution to the gravitational interaction is of the Yukawa type, if oc GMe~"^^ /rc^, still not 
of the MONDian form, notably because its global factor is M instead of a/M. [Note that 
the Yukawa force is proportional to dr^p oc GMe~"^^{l/r'^ + m/r), and therefore includes a 
MOND-like 1/r contribution. However, it dominates the main 1/r^ term only for mr ^ 1, 
i.e., precisely when the exponential e~"^^ makes the whole force quickly tend towards zero.] 
A "quintessence" -hke potential, whose minimum occurs at — oo, seems better because it 
allows us to build a logarithmic gravitational potential. Indeed, if V{(p) = —2a'^e~'"^, where 
a and b are two constants, then (f = (2/6) ln(a6r) is a solution of the vacuum field equation 
Aip = V'{(p). However, not only this potential is unbounded by below, thereby spoiling the 
stability of the model, but if is also multiplied by a constant 2/6 (independent of the matter 
source) instead of being proportional to a/M. 



A generalization of f{R) theories has been considered in Refs. [8J, |95|, where the 
Lagrangian density is given by some function of the scalar curvature and its iterated 
d'Alembertian, f{R, OR, . . . , □"/?). The theory can be shown to be generically equivalent 
to a tensor- multi-scalar theory |91i] involving n + 1 scalar fields. The signs of their respective 
kinetic terms depends on the numerical constants entering the function f{R, OR, . . . , □"/?), 
and there is thus no guarantee that they all carry positive energy in the general case. How- 
ever, well chosen functions / do define stable tensor-multi-scalar theories. 

Although it seems difficult to reproduce MOND dynamics with tensor-mono-scalar mod- 
els, general tensor-multi-scalar theories might actually succeed. Indeed a particular tensor- 
bi-scalar model, called "phase coupling gravity" (PCG), has been shown to reproduce the 
TuUy-Fisher law in some ap propriate regime (there exists a critical mass above which MOND 
phenomenology arises) [96|, l97|. However this model is marginally ruled out experimentally 
and moreover needs a potential which is unbounded by below, thereby spoiling the 



stability of the model. This kind of theory has then been promoted to a "tensor-vector-bi- 



scalar" model [33|, that we shall examine in more detail in Sec. ITTO 

Among stable scalar-tensor theories, let us finally mention those involving a coupling of 
the scalar field to the Gauss-Bonnet topological invariant, i.e., involving a term of the form 
f{^p) X GB in the gravitational action. [Note that the function of is a priori free, but 
that the Gauss-Bonnet term must appear linearly.] Contrary to the higher-order theories 
discussed in Sec. Ill CI above, such models do not always cont ain g host degr ees of freedom. 
They have been studied in various contexts (see e.g. (isl, [9^, 100 . 101 . 102l |). but not with 



the aim of reproducing the MOND phenomenology. Our own investigations showed that it 
is a priori promising, because the Riemann tensor R^upa generated by a body of baryonic 
mass M does not vanish outside it, and therefore gives us a local access to M. However, a 
second information is also needed to separate M from the radial dependence, and we did 
not find any simple way to generate a potential ^/GMao Inr. [It is actually possible if one 
defines a nonstandard kinetic term for the scalar field, as in Sec. Ill El below, but the scalar- 
Gauss-Bonnet coupling is not necessary in such a case.] A related but much simpler idea is 
explored in Sec. |V] below. 
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E. K-essence or RAQUAL models 



The above scalar-tensor models may involve three functions of the scalar field, namely 
the potential V{ip), the matter-scalar couphng function A{{p), and the possible scalar-Gauss- 
Bonnet coupling /(v?) X GB. Their generalization to n scalar fields j9l| may also in- 
volve a n X n symmetric matrix 'jabi'P'^) depending on them and defining their kinetic term 
g'i'^lab{'^^)diiV"'du^^ ■ Their phenomenology is however similar to the single scalar case, at 
least if all of them carry positive energy and their potential is bounded by below, to ensure 
the stability of the theory [91i|. On the other hand, a significantly different physics arises 
if one considers more general kinetic terms of the form f{s,ip), where s = g^^d^ipdyip is 
the standard one. [Beware that the literature often uses the notation X = s/2, which ac- 
tually does not change the following discussion.] Such models, without any matter-scalar 
coupling, have been studied in the cosmological context under the name of "k-infiation" [24] 



or "k-essence" [251, |26|, the letter k meaning that their dynamics is kinetic dominated, as 
opposed to quintessence models in which the potential V{ip) plays a crucial role. When a 
matter-scalar coupling A{ip) is also assumed, such models have been studied to reproduce the 
MOND dynamics under the name of "Relativistic AQUAdratic Lagrangians" (RAQUAL) 
27 . Their action is thus of the form 



^=4^/ - ^f^'^^) - ^^^)] + ^rn^Uerii^rg^u = A\^)g;,], (2.13) 

with s = gi^'^d^ifd^ip, and they can be considered as natural generalizations of the simplest 
f{R) model, Eq. (12.121) . Note that the potential V{!f) may be reabsorbed within the general 
function f{s,!f), and that it is therefore unnecessary. However, since one often considers 
functions f{s) of the kinetic term alone, it remains convenient to keep an explicit potential 
V{(f). Reference 13] proved that an appropriate nonlinear function f{s,ip) allows us to 
reproduce the MOND gravitational potential ^/GMaolnr, but that some difficulties remain, 
as discussed below and in Sees. IIIFI and IIIII 

Several conditions must be imposed on this function f{s, ip) to guarantee the consistency 
of the field theory. For any real value of s G [—00, +00], one must have 

(a) f(s,<^) >0, 

(b) 2sf"is,^) + f'is,^)>0, 

where a prime denotes derivation with respect to s, i.e., /' = df{s,(p)/ds. The condition 
f'{s) > 0, implied by (a), is necessary for the Hamiltonian to be bounded by below. Indeed, 
up to a global factor c^/SirG that we do not write to simplify the discussion, the contribution 
of the scalar field to this Hamiltonian reads H = 2{doip)'^ f (s) + f{s) in a locally inertial 
frame. If there existed a value s_ such that /'(s_) < 0, then H could be made arbitrarily 
large and negative by choosing diverging values of {doipY and {diipY such that — (SqV?)^ + 
{diipY = s_. On the other hand, conditions (a) and (b) are necessary and sufficient for the 
field equation to be always hyperbolic, so that the Cauchy problem is well posed for the 
scalar field. Indeed, it reads 
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where G^^ = fg^'^ + ^fV^fV^ip plays the role of an effective metric in which (p propagates. 
The lowest eigenvalue of G'^'^ is negative (i.e., defining a consistent time) only if condition 
(b) is satisfied, and it is easy to check^ that the three others are then positive (defining 
spatial dimensions) if (a) is also satisfied.^ As a quick check of condition (b), one may for 
instance consider the particular case of a homogeneous scalar field {di^p = 0) in a locally 
inertial frame (5?*^ = t]^^): The fact that we need (7°° = — /' + 2f" x (— s) < then 
immediately yields inequality (b). Finally, conditions (a) and (b) together suffice for the 
Hamiltonian to be bounded by below, provided the function f{s,(p) is analytic and f{s = 
0, if) is itself bounded by below for any p. Indeed, still in a locally inertial frame, we know 
that s = —{doip)'^ + (diip)'^ > —{doip)'^. Together with condition (a), we can thus conclude 
that the Hamiltonian H = 2((9oV5)^/' + / is greater than —2s/' + /, which is known to be a 
decreasing function of s because of condition (b). Therefore, for any s < 0, the Hamiltonian 
is necessarily greater than f{s = 0,ip). On the other hand, for s > 0, we also know that 
H > f > f{s = 0,{p), first because (dop)"^ > and then because condition (a) implies that 
/ is an increasing function of s. Although we used the hyperbolicity condition (b) to derive 
that the Hamiltonian is bounded by below, note that condition (b) is not implied by this 
boundedness, on the contrary. Nevertheless, it reappears when considering the propagation 



velocity of perturbations around a background [2J], = /'/ {2s f" + /'). Such perturbations 
are unstable if < 0, and one may thus think that condition (b) should be implied by 
the boundedness by below of the Hamiltonian. However, the energy of the perturbations is 
only part of the total Hamiltonian, which also includes the energies of the background and 
its interactions with the perturbations. A total Hamiltonian bounded by below therefore 
does not suffice to guarantee the stability of perturbations; the hyperbolicity of the field 
equations, condition (b), is al so crucial. 
In Ref. [l03 (see also [lo|), 

a third condition was underlined for such k-essence models, 

besides (a) and (b) above: 
(c) f'is) < 0. 

This extra condition is necessary for the causal cone of the scalar field to remain inside 
the light cone defined by the spacetime metric g^^, i.e., to avoid superluminal propagation. 
The simplest way to derive this inequality is to consider a null vector k'^ with respect 
to the metric G^^ = il/f')[g*^^ - 2f"d^^d^ip/{2sf" + /')], defined as the inverse of the 
effective metric G^^ entering Eq. fl2.14p above. [We assume here that /' never vanishes.] 
The equation G^^kf'k" = then gives g*^^k^'k'' = 2f"{k>'df,pf/{2sf" + /'), and obviously 
implies g*^^k^k'^ < if conditions (b) and (c) are satisfied. The causal cone of the scalar field 
is thus timelike or null with respect to the Einstein metric g*^^^. Such a condition (c) would 
be extremely constraining, since it would rule out notably any monomial f{s) = s", which 
would violate either (a) or (c) depending on the parity of n. However, let us underline that 



^ The simplest calculation consist in diagonalizing the matrix C^^g*^, instead of the tensor G^'^, and impose 
that its four eigenvalues are positive. This also ensures that there exists a coordinate system in which the 
time directions defined by G^i/ and 5*^ coincide. 

^ One could also consider the degenerate case where /'(s) vanishes on some hypersurfaces. On them, the 
effective metric can then be put in the form diag(— 1, 0, 0, 0) in some appropriate basis, and thus define 
a consistent time, if and only if /"(s) < 0. Strictly speaking, the scalar field equation is not hyperbolic 
in that case but it has nevertheless a well-posed Cauchy problem. Using condition (b), we can conclude 
that /'(s) may vanish only at negative values of s. 
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causal cones 




FIG. 1: In a theory where different fields have different causal cones, it suffices that their union 
be embedded in a wider cone for local causality to be satisfied. Initial data for all the fields 
simultaneously may be specified on a surface exterior to the wider cone, i.e., spacelike with respect 
to each cone. If the topology of spacetime is such that there does not exist any CTC with respect 
to the wider cone, then causality is preserved, although some fields may propagate faster than light 
(i.e., faster than electromagnetic waves, a mere particular case of matter field). 



(c) is actually not required for the consistency of the model. Indeed, even if propagates 
superluminally, this does not ruin the causality of the model. Thanks to condition (b), the 
causal cone of never opens totally, i.e., always admits a class of Cauchy surfaces which 
can be defined in its exterior. In the general situation where different fields have different 
causal cones, it suffices that their union still admits a nonvanishing exterior where one may 
consistently define initial data for all of them simultaneously; see Fig. [TJ This becomes quite 
clear by shifting our viewpoint. Let us assume that the union of all causal cones may be 
embedded in the interior of a wider cone which never totally opens, and let us define space 
and time with respect to this exterior cone. Since all fields propagate within this cone, there 
is no more causal problem than in special relativity where particle worldlines lie within the 
standard lightcone. Causal pathologies are still possible, like in GR itself which admits 
solutions containing closed timelike curves (CTC), like the Godel universe or the interior of 
the Kerr solution.^ Since the global topology of the Universe is not imposed by our local field 
equations, it is necessary to assume that it does not involve any CTC to ensure causality. In 
our example of different fields having different causal cones, one thus just needs to assume 
that spacetime does not admit any CTC with respect to the wider cone surrounding all of 
them. Caus al pa radoxes due to superluminal propagation, like the nice one exhibited in 



Sec. III. 2 of |104l |. are thus more a mat ter of global assumptions than a problem with local 



field equations.^ We refer to Ref. jl06l | for a more detailed discussion of this subtle issue 



^ It has even been proven that CTCs can form in G R fro m smooth initial data, while assuming asymptotic 

spatial flatne ss and various energy conditions; see |105l | and references therein. 
^ In Ref. 10J|, a CTC is exhibited by considering a background involving two bubbles of a scalar field 



with a fast relative motion. If one assumes that all clocks can be globally synchronized by using only 
elect romagnetic waves (for instance that spacetime is Minkowskian with respect to 17*^ = 77^^, like in 



lOJl), then such a CTC does exist. However, although this hypothesis sounds extremely natural, there 



also exist spacetimes in which there are no CTCs for any field, and in such a case, light within a given 
bubble of scalar field does not define the same synchronization as within the other bubble. The CTC 
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(see also Sec. IIII CI below and the very recent article 107l |). Reference jl08 | claims that 



microscopically Lorentz-invariant particles cannot give rise to superluminal signals, but this 
conclusion does not take into account their possible self-interactions via their kinetic term, 
like f{s) in our present case of k-essence models fl2.13p . As shown above, superluminal 
signals do occur when f"{s) > 0, although the theory is microscopically Lorentz-invariant 
— and causal when the hyperbolicity conditions (a) and (b) are satisfied. The causal cone 
of the scalar field (p is either interior or exterior to the light cone defined by (7*^ (depending 
on the local value of s), but it always remains a cone thanks to these conditions. At each 
spacetime point, any surface exterior to the wider of these two cones may be used as a 
Cauchy surface to impose initial data. Of course, a surface lying between the two cones 
behaves as a spacelike one for the thinner cone but as a ti melik e one for the thicker. Some 
causal paradoxes discussed in the literature (including in jl03l |) are actually based on the 
improper use of such intermediate surfaces. They are clearly not consistent Cauchy surfaces 
for the fastest field (i.e., the wider cone), so that data are necessarily constrained on them. 
The conclusion of such paradoxes, namely that initial data are constrained, is thus actually 
hidden in the use of these intermediate surfaces. 



F. MOND as a RAQUAL model 



As mentioned above, Ref. [27[ proved that RAQUAL models 02.131) may reproduce the 
MOND gravitational potential. Indeed, if /(s) is assumed not to depend explicitly on (p, if 
the potential V{(p) vanishes, and if one chooses an exponential (Brans-Dicke-like) matter- 
scalar coupling function A{ip) = exp{aip), Eq. fl2.14p reduces to 

V;[/(.)Vi:v.] = -^ar*, (2.15) 

where T* = g*^JTi^'^ denotes the trace of the matter energy-momentum tensor T^'^ = 
(2c/ V— (7*) 5 Scatter 1 5 9%- In order to recover the right Newtonian limit for large accelera- 
tions a > ao, one may follow Ref. ^27j and impose^ f'{s) 1 for large positive values of s. In 
such a case, the scalar field reads cp ~ —aGM/rc^ near a body of mass M, and test particles 
feel an extra potential apc^ ~ —a^GM/r adding up to the standard gravitational potential 
—GM/r mediated by the spin-2 interaction (i.e., by the Einstein metric g*^j^y)- The total 
potential — G(l -|- a^)M/r is thus of the Newtonian t ype , with a renormalized effective grav- 



itational constant Gefi = G{1 + a^) (see for instance |9l| and line (a) of Fig. [5] in Sec. Illllbe- 
low). Precision tests in the solar system are however sensitive to post-Newtonian corrections, 
and they prove that the scalar contribution must be negligible. Indeed, the parametrized 
post-Newtonian (PPN) parameter 7^™ assumes the value 7^™ = 1 — 2a^/(l + a^) in the 
present (conformally-coupled) scalar-tensor framework [see Refs. [H, [9l| and Eq. (03]) be- 
low], and the impressive experimental bound [7^^^ — 1| < 2 x 10~^ obtained in [90] from 



104l | would then become an artifact of gluing independent coordinate systems at the edges 



constructed in 
of the bubbles. 

® Note that imposing f'{s) const, is strictly equivalent, since the arbitrary constant can be absorbed in 
a redefinition of the matter-scalar coupling constant a. In the following, the experimental constraints on 
a that we quote correspond to the choice /'(s) ^ 1 in the Newtonian regime. 
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FIG. 2: A simple function /'(s) reproducing the MOND dynamics for small s (i.e., large distances) 
and the Newtonian one for large s (i.e., small distances). 



the observation of the Cassini spacecraft implies therefore 

< 10-\ (2.16) 

Instead of imposing f'{s) — ^ 1 for large values of s, one may also recover the right Newtonian 
limit by choosing a decreasing function f'{s), such that the scalar contribution is even 
smaller. Condition (b) should however still be satisfied, thereby constraining the slope of 
fis). 

But the crucial feature of Eq. (12.151) is that it also allows us to reproduce the MOND 
potential for small accelerations a < qq. Indeed, if /'(s) ~ io\/s for small and positive 
values of s, where io is a constant having the dimension of a length, one gets {dr^Y ~ 
aGM I i^c^r^ near a spherical body of mass M. Test particles therefore feel, in addition 
to the usual Newtonian potential —GM/r, an extra potential a^pc^ ~ a/ ct'GMcp- j In r 
which reduces to the MOND one ^/GMa^ Inr for = ot'(? ja^. A simple way to connect 
this MOND limit to the above Newtonian one is to choose for instance /'(s) = a/I/a/I + s, 
where s = a^c^s/a^ is dimensionless. The shape of this function is illustrated in Fig. [21 
By integrating this expression, one can thus conclude that the precise function f{s) = 

(a^/aV) V^(l + s) -sinh-^(v^) , or any other one having the same asymptotic behaviors 

for large and small positive values of s, allows us to reproduce the MOND dynamics. 

However, the authors of Ref. [23| noticed that the scalar field propagates superluminally 
in the MOND (small positive s) regime, since f"{s) ~ \io/\/s > contradicts condition 
(c) above. This is the reason why they discarded this model, although we underlined above 
that such a superluminal propagation actually does not threaten causality, provided the 
hyperbolicity conditions (a) and (b) are satisfied. In fact, the only experimental constraints 
we have on the existence of several causal cones is in favor of superluminal fields! As we 
will discuss in Sec. IIV C\ if light (and thus matter) travelled faster than some field to which 
it couples, then it would emit Cherenkov radiation of this field. High-energy cosmic rays 
would thus be significantly suppressed. 

It thus seems that the above RAQUAL model reproducing the MOND dynamics is a 
priori consistent. However, it presents several other difficulties. The main one has been 
immediately recognized and addressed in the literature: Because of the specific form of the 
matter-scalar coupling assumed in action (12.131) . this model does not reproduce the observed 
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light deflection by "dark matter". We will devote Sec. UTIl below to this crucial problem. Let 
us here mention the other problems that we noticed, and which do not seem to have been 
discussed in the literature. 

First of all, the above function f{s) is clearly not deflned for negative values of s. One 
may try to replace s by its absolute value \s\, and also multiply globally /(s) by the sign of 
s in order to still satisfy conditions (a) and (b). However, this would not cure the serious 
problem occurring at s = 0. Indeed, since f{s) (2a^c^/3ao)sA/|s| for small values of 
s, the strict inequality (b) is not satisfled at s = 0. In other words, the scalar fleld equa- 
tion fl2.14p or fl2.15p is no longer hyperbolic on the hypersurfaces where s vanishes, and the 
scalar degree of freedom cannot cross them. Since s ~ {di^Y — "when considering the local 
physics of clustered matter, but s ~ —0^ < when considering the cosmological evolution 
of the Universe, there always exist such singular surfaces around clusters. Therefore, this 
model cannot describe consistently both cosmology and galaxies, unless independent solu- 
tions are glued by hand on the singular hypersurfaces s = 0. However, a simple cure to this 
discontinuity would be to consider, for small values of s, a function 



f\s)=e+^/\s\, (2.17) 

where e is a small dimensionless positive number. Then both conditions (a) and (b) would 
obviously be satisfled even at s = 0. In other words, the RAQUAL model (I2.13p with 

V{(p) = 0, A{{p) = exp{a(p) and 

fis) = es + [VmrW) - sinh-^ (v^)] , (2.18) 

where s = a^c'^s/aQ, does reproduce the MOND dynamics and is free of mathematical 
inconsistencies. Of course, there is no reason why this speciflc function, for negative values 
of s, should reproduce the right cosmological behavior. One may obviously look for other 
functions of s < connecting smoothly to fl2.18p for small values of \s\, provided conditions 
(a) and (b) above remain always satisfled. Therefore, Eq. (12.180 should just be considered 
as an example of a mathematically consistent RAQUAL model. 

A second difficulty is that the above model is rather fine tuned, since it needs the intro- 
duction of a small dimensionless constant e besides Milgrom's MOND acceleration constant 
ao ~ 1.2 X 10~^°m.s~^. Indeed, the presence of e in Eqs. fl2.17l) . fl2.18p notably implies 
that Newtonian gravity is recovered at very large distances, the MOND regime manifesting 



only at intermediate ranges [106|. When s — > 0, the derivative f'{s) entering Eq. fl2.15p 
tends to e, so that the scalar fleld reads ip ~ —{a/e)GM/rc^ faraway from a body of 
mass M. In this regime, the total gravitational potential felt by a test particle reads 
—G{1 + a'^/e)M/r, and is therefore of the Newtonian form with a renormalized gravita- 
tional constant Goo = ^(1 + a'^/e), where the subscript refers to the fact that the above 
form of the potential holds exactly when r — cx) (remind that, if r — > 0, the effective grav- 
itational constant is given by Gefr = G(l + a^)). Let us compute the range of distances 
r for which the MOND force dominates the Newtonian contribution. One needs of course 
f'{s) oc ^/s and thus e <^ ^/I. On the other hand, the MOND force dominates the Newto- 
nian one if r > a/GM/oo. Using s ^ {dripf ^ GMao/ a^c^r'^ , we thus flnd that the MOND 
force dominates within the following range of distances ^J CM/ ao -C r <^ {a^ / e) ^/gWJoq. 
Since solar system tests impose < 10 ~^, Eq. (12.160 . and since rotation curves of galaxies 
may be flat up to r ~ lO^GM/oo flOgj . one needs therefore e -C 10 ^. This illustrates the 



fine tuning required to define a consistent RAQUAL model even for s — 0. 
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On the other hand, Eq. (12.171) has two great advantages, besides the fact that it cures 
the singularity at s = 0. First, the fact that the field assumes a Newtonian form at 
large distances, r > {a'^ /e)^^/GMja(), ensures that a body of mass M does not create a 
static (scalar) field with infinite energy, contrary to the case where f'{s) ~ y/s for small 
23]. Secondly, the presence of e naturally leads to a theory where the local constant 



of gravity Gefr = G{1 + a^) may differ from Goo, i-e., the one in regions of spacetime 
where the scalar field gradient is small. This may have some interesting applications in 
cosmology. Indeed, provided that \/\^ ^ e, i.e., when 0^ < eal/a^c^ in a cosmologi- 
cal context, a homogeneous Universe is described by Friedman equations^ but the bary- 
onic matter contributes as GooPbaryons- The effective density of matter reads therefore 
PbaryonsGoo/Geff = Pbaryons(l + a'^ / s) / {1 + Qf^) f» Pbaryons(l + a'^/e). In other words, bary- 
onic matter may contribute effectively as if there were true (exotic) dark matter in that 
regime, provided that we tune the free parameter e so that a'^/e ~ 10, which is consistent 
with galactic dynamics. Since the cosmology of such models goes beyond the scope of the 
present paper, we will not investigate further this phenomenology. Let us however stress 
that there is here a potential threat for MOND-like theories. Indeed, bounds may be found 
on the ratio Goo/Ges from the analysis of primordial nucleosynthesis, since the speed of the 
expansion of the Universe, and therefore the relative proportion of hydrogen and helium in 
the present Uni verse, dep ends on this ratio. W ithin the framework of the so-called Einstein- 
aether theories lid . Ill ], it was shown in Ref. 112] that one must have \Goo/Gcs — 1] ~ 1/8 



to be consistent with the data. Of course, in the above RAQUAL model of MOND, such 
a constraint only exists if the time derivative of the scalar field is small at the epoch of 
nucleosynthesis. In that case, the above bound would read a'^/e < 1/8, and the MOND 
regime would therefore not exist at all (there is no more room for the flattening of rotation 
curves, as shown in the above paragraph). A solution would of course to impose a very small 
value of e such that ip^ ^ eal/a^c^ at the nucleosynthesis epoch. However one should study 
the cosmology of such models to investigate in more detail this effect and make quantitative 
estimates; we leave this point for future work.^° 

Finally, the third (and the most serious) problem of such RAQUAL models of MOND, 
besides the one of light deflection discussed in Sec. [1111 below, is that the MOND force 
oc y/GMao/r starts manifesting at quite small distances. Indeed, the MOND potential 
oc V GMao In r appears at values of field gradient s such that f'{s) ~ v^. This typically 
happens at s ^ 1 for natural functions /' interpolating between the MOND and the Newto- 
nian regimes, see Fig. [2l Using s ^ 1 we find that the MOND potential appears at distance 
r ^ rtrans f^om a body of mass M, where 



9 GM 

rtrans = • (2.19) 

V O'O 

Since, in the present (conformally-coupled) scalar-tensor framework, the matter-scalar cou- 
pling constant a is severely constrained by solar-system tests, Eq. (12.161) . we conclude that 
the extra MOND potential starts manifesting at a radius rtrans = (y^ \/ GM^j a^) < 0.1 AU 



^ Beware that the density and the pressure of the scalar field must also be taken into account. 

° Note that this phenomenon is very likely to occur also in extended models we consider below (TeVeS-like 

theories), since they notably involve a k-essence scalar field whose kinetic term also needs to be cured at 

s = 0; see Sec. HITCl 
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FIG. 3: Fine-tuned function f'{s) such that Newtonian and post-Newtonian predictions are not 
spoiled in the solar system, although the MOND dynamics is predicted at large distances. The 
right panel displays the quite unnatural contribution of the scalar field to the acceleration of a test 
mass, as a function of its distance with respect to the Sun. 



in the solar system, where Mq denotes the mass of the Sun. In other words, all plan- 
ets (including Mercury, at 0.4 AU from the Sun) should feel a total gravitational potential 
—GMQ/r + y/GMQao Inr, the first term coming from the spin-2 field g^^^, and the second one 
from the scalar field which is already in its MOND regime. This extra potential y^GMQao In r 
(or force in 1/r) leads to deviations f rom the Newtonian behavior and is tightly constrained 
by tests of Kepler's third law; see jll3 | and Sec. IV CI below. Although this MONDian 
anomalous force remains numerically small with respect to the Newtonian one within the 
solar system, it is large with respect to post-Newtonian corrections oc and its exis- 

tence is definitely ruled out experimentally. In conclusion, although the RAQUAL model 
of Ref. [23] or its refinement (12.181) do predict a Newtonian behavior at distances r — 
(i.e., s — >■ oo), planets are far enough to be already in the MOND regime! Note that this 
surprizing problem comes from the experimental bound (I2.16p . telling us that the scalar 
field must be very weakly coupled to matter to reproduce GR within the solar system. If 
one did not take it into account, by setting a = 1 for instance, one would easily reproduce 
Newtonian gravity up to a/ GMq / ~ 7000 AU, and the MOND dynamics beyond. A 
reasonable value e ~ 1/10 would also suffice for the refined model (I2.18p . However, this 
theory would be inconsistent by 5 orders of magnitude with post-Newtonian tests in the 
solar system, because the scalar field would be much too strongly coupled to matter. 

A possible solution to the above problem would be to fine tune the function f{s) even 
further. In order to get the MOND regime for r > a/GM/oo, as required by galaxy rotation 
curves, one would need f'{s) ~ y/i for s ^ < 10~^°. On the other hand, in order to 
obtain the Newtonian regime within the solar system, say for r ^ rmax ~ 20 or 30 AU, 
one would need f'{s) ~ 1 for s ^ (a^GM0/aor^ax)^ ~ 10"^° (this second occurrence of 
10~^° is a numerical coincidence). Therefore, there would exist a brutal transition between 
the MOND and Newtonian regimes around s ~ 10"^''. Not only the introduction of such 
a small dimensionless number would be quite unnatural, but this model would also predict 
that the anomalous acceleration caused by the scalar field remains approximately equal to 
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FIG. 4: Typical shape of a function f'{s) such that scalar-field deviations from GR are oc r", n > 0, 
in the large-s limit. 



the constant^^ ao between 30 and 7000 AU. As illustrated in Fig. [3l this would be a way 
to reconcile the MOND acceleration y/GMa^/r at large distances with the experimentally 
small Newtonian contribution a^GM/r"^ of the scalar field at small distances. Although 
this is not yet excluded experimentally, it would however suffice to improve by one order 
of magnitude the post-Newtonian constraint on to rule out such a fine-tuned model 



(the planned astrometric expe riment G AIA |114| | should reach the 10 level for a , and 



the proposed LATOR mission fTT5l.[TT6j should even reach the 10 ^ level). Therefore, one 
should not consider it seriously. 

An a priori better solution to the above problem would be to recall that the Newtonian 
limit /'(s) ^ 1 is actually unnecessary. Since the metric g*^^^ already generates a Newtonian 
potential —GM/r, it suffices that the scalar field contribution remain negligible (even at the 
post-Newtonian level) in the solar system. One may try a function f'{s) whose shape looks 
like the one displayed in Fig. HJ for instance 



f'{s)=e + ^ ' ' , (2.20) 

(l + |s|)'+'/" V ) 

where n is a positive constant and as before s = a^c^s/oQ. Then the transition occurs again 
at s ~ 1, i.e., around rtrans = a^^GM/ao, but Eq. (12.151) shows that the force mediated by 
the scalar field reads adrfc^ = (ao/«^) / ''^trws)^ in the large-s limit [assuming that the 
e of Eq. (I2.20p has a negligible infiuence]. Therefore, even if a = 1, so that the Newton- 
MOND transition occur at rtrans = \/ GM/ oq as expected, the anomalous scalar force would 
be negligible with respect to post-Newtonian relativistic effects for r -C rtrans and n large 
enough. However, the above scalar field drf oc r" happens not to be a solution of Eqs. ( 12.15^ 
and (I2.20p for r < rtrans, where e must actually dominate. In such a case, the scalar force 
takes the Newtonian form adr^pc^ = a^GM/r'^e and is even increased by a factor 1/e with 
respect to model (I2.17p . This suffices to rule out the class of models (12.201) . One may also 



Note that the MOND constant ao is too small by a factor 7 to account for the anomalous acceleration 
of the two Pioneer spacecrafts. Actually, we will also see in Sec. IV CI below that there exists a crucial 
difference between the MOND dynamics and the Pioneer anomaly. 
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notice that for e ^ 0, they do not satisfy the hyperbohcity condition (b) in the large-s hmit 
unless n = oo; but even for this hmiting case n = oo, the scalar force is not negligible at 
small distances. 

In conclusion, the above discussion illustrates that RAQUAL models are severely con- 
strained, although they involve a free function f{s) defining the kinetic term of the scalar 
field. Contrary to some fears in the literature, the possible superluminal propagations do 
not threaten causality, and the two conditions (a) and (b) are the only ones which must 
be imposed to guarantee the field theory's consistency. For instance, monomials f{s) = s"' 
are allowed if n is positive and odd [except on the possible hypersurfaces where s vanishes, 
which would violate the strict inequality (b)]. However, when taking into account simulta- 
neously these two consistency conditions and experimental constraints, it seems difficult to 
reproduce the MOND dynamics at large distances without spoiling the Newtonian and post- 
Newtonian limits in the solar system. It seems necessary to consider unnatural functions 
f{s) involving small dimensionless numbers. 

In the next section, we will address the problem of light deflection, and recall the solution 
which has been devised in the literature. This solution will at the same time release the 
experimental constraint on the matter-scalar coupling constant a obtained from accurate 
measurement of the post-Newtonian parameter 7^^^ within the solar system, and therefore 
cure the above fine-tuning problem. On the other hand, we will show in Sec. IIVDI that 
the analysis of binary-pulsar data also imposes a bound on a, so that the problem of the 
fine-tuning of the free function / will anyway remain. 



III. THE PROBLEM OF LIGHT DEFLECTION 
A. Conformally-coupled scalar field 

Both scalar-tensor theories fl2.12p and RAQUAL models (12.131) assume a conformal re- 
lation Qfj^i, = A'^[ip)g^^ between the physical and Einstein metrics. The inverse metrics are 
thus related by g'^" = A~'^{(p)gl^'^ , and their determinants read g = A'^'^{ip)g^: in d spacetime 
dimensions. Since all matter fields, including gauge bosons, are assumed to be minimally 
coupled to (jfj,!^, this leads to a simple prediction for the behavior of light in such models. 
Indeed, the action of electromagnetism is conformal invariant in d = 4 dimensions: 

= / ^ [A-\^)9r] F,.F,^ 

= / V ^9^'97F,.F,.- (3.1) 

Therefore, light is only coupled to the spin-2 field (7*^^, but does not feel at all the presence 
of the scalar field (p. In terms of Feynman diagrams, there exist nonvanishing vertices 
connecting one or several gravitons to two photon lines, but the similar vertices connecting 
photons to scalar lines all vanish; this is illustrated in line (b) of Fig. [5l It is then obvious 
that light behaves strictly as in GR, in a geometry described by the Einstein metric f/*^. In 
particular, the light deflection angle caused by a spherical body of mass M must be given 
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FIG. 5: Feynman diagrams in scalar-tensor theories, where straight, curly and wavy lines represent 
respectively the scalar field, gravitons and photons. Matter sources are represented by blobs, 
(a) Diagrammatic interpretation of the effective gravitational constant Ges = G'(l-l-ao), where each 
vertex connecting matter to one scalar line involves a factor oq. (b) Photons are directly coupled 
to gravitons but not to the scalar field, (c) Photons feel nevertheless the scalar field indirectly, via 
its infiuence on gravitons: The energy-momentum tensor of the scalar field generates a curvature 
of the Einstein metric g* in which electromagnetic waves propagate. 



by the same expression as in GR (at lowest order) 

M = (3^2) 

where h denotes the impact parameter of the light ray. An even simpler way to prove that 
light propagates in the Einstein metric (7*^, without feeling the scalar field, is to note that 
its geodesic equation reads g^^dx^dx^ = in the eikonal approximation. Dividing by the 
nonvanishing factor this equation implies g*^^dx^dx^ = 0, giving thus the standard 

geodesic equation for null rays in the metric g*^^. [It is interesting to note that this second 
reasoning would remain valid even if the action of electromagnetism were not conformal 
invariant, for instance in dimension d ^ 4, or even if one multiplied it by an explicit scalar- 
photon coupling function B'^{ip), thereby violating the weak equivalence principle. Then it is 
straightforward to prove that the nonvanishing scalar-photon vertices would only affect the 
amplitude of electromagnetic waves, in the eikonal approximation, but not their polarization 
nor their trajectory. One of us (G.E.F.) disc ussed this result with B. Bertotti several years 



ago, but it has been recently rediscovered in [117 



However, the crucial difference with GR is that massive matter does feel the scalar field, 
via the matter-scalar coupling function A{ip). One may expand this function around the 
background value (po of the scalar field as 

In A{ip) = const. + ao(v2 - V'o) + ^Poiv - V'o)^ H , (3.3) 

where ao, Po, are dimensionless constants. In usual scalar-tensor theories (12.121) . as 
well as in the Newtonian regime of RAQUAL models (I2.13p . the lowest-order contribution 
of the scalar field to the gravitational potential felt by a test mass reads then —a^GM/r 
[l|, 9l[. This generalizes the results recalled in Sec. Ill Fl above in the particular case of an 



exponential coupling function A{ip) = e°"^. This spin-0 contribution adds up to the standard 
Newtonian potential —GM/r caused by the spin-2 interaction (i.e., via the Einstein metric 
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(7^^), so that the total gravitational potential remains of the Newtonian form but with an 
effective gravitational constant Gcs = G{1 + ao)', see line (a) of Fig. [5l In other words, when 
one determines the mass of the Sun by analyzing the orbits of the planets, one actually 
measures GesM instead of the bare GM. Therefore, if one expresses the light deflection 
angle (I3.2p in terms of the observable quantity GesM, one gets 

- {l + al)bc^ - -b^- ^^-^^ 

For a given value of the observed constant GesM, general relativity would have predicted 
4:GesM/bc^. In conclusion, the light deflection angle is actually smaller in scalar-tensor 
theories than in GR, although photons are strictly decoupled from the scalar field! This 
paradoxical result is simply due to the fact that we need to compare the measured angle 
with another observable, G^ffM, which is influenced by the scalar field. 

This reduction of light deflection in scalar-tensor theories has been immediately recog- 



nized as a serious problem for MOND-like field theories [28|, |29l]. Indeed, if one interprets 



data within GR, the presence of dark m atter sug gested by galaxy rotation curves is con- 
firmed by weak lensing observations [ iTsl . liol . 120l | . Therefore, in a modified-gravity theory 
avoiding the assumption of dark matter, one should predict simultaneously a larger gravita- 
tional potential and a larger light deflection. However, the interpretation of data is slightly 
more subtle than in the solar system, because there exist at least three different notions of 
mass which may be defined for a galaxy: (i) its baryonic mass Mb, assumed to be propor- 
tional to its luminous mass Ml; (ii) its total mass (baryonic plus "dark") evaluated from 
its rotation curve, that we will denote M^°^ ; and (iii) its total mass evaluated from weak 
lensing, that we will denote Mll^. The only conclusion driven from the above result fl3.4p 
is that the deflecting mass is smaller than the gravitational one in scalar-tensor theories: 

MlT < Mi:l. (3.5) 

Note that this conclusion is also valid in the class of RAQUAL models reproducing the 
MOND dynamics that we discussed in Sec. Ill Fl above. Indeed, the MOND force was caused 
by the exchange of a scalar particle [cf. line (a) of Fig. [5], i.e., by the conformal factor A^{ip) 
relating the physical metric g^j,^ to the Einstein one g*^^. Since light only feels the latter, 
the total mass evaluated from weak lensing is obviously smaller that the one evaluated from 
rotation curves. 

On the other hand, it should be underlined that Eq. fl3.4p does not imply M\^^ < Mb, 
contrary to some erroneous or imprecise claims in the literature (see, e.g., [1^). Indeed, 
the stress-energy tensor of the scalar field does contribute to the curvature of the Einstein 
metric ^f*^, so that light deflection may be larger than if only baryonic matter was present. 
This is illustrated in terms of Feynman diagrams in line (c) of Fig. [5j Although the scalar- 
photon vertex vanishes, light may nevertheless feel the scalar field indirectly via a graviton 
exchange. The reason why this does not appear in Eqs. (13. 2p or (13. 4p above is because we 
derived them at lowest order in powers of and in usual scalar-tensor theories or in the 
Newtonian regime of RAQUAL models. But when considering their MOND regime, one 
cannot any longer perform a naive post-Newtonian expansion. Indeed, the stress-energy of 
the scalar field may become non-negligible with respect to that of baryonic matter, and even 
dominate it, although this scalar field has been generated by its coupling to matter. This 
mere sentence gives us a hint that such a model may be unstable, but one should note that 
Mix > is indeed possible in RAQUAL models. 
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We will provide an explicit example in Sec. IIIIBI below, although it does not reproduce 
the right MOND dynamics. Its aim is just to underline that the derivation of the inequality 
^tot'^ — proposed in Ref. j29| is incorrect. It is instructive to locate where its reasoning 
fails. It writes the deflection angle in terms of the component of the total stress-energy 
tensor, and decomposes it in two contributions, due to matter and to the scalar field. It 
happens that the latter is negative. However, what was identified as the matter contribution 
to happens not to vanish outside matter. Therefore, the decomposition was flawed, and 
part of the positive contribution of the scalar field to light deflection was actually attributed 
to matter. The simplest way to understand that a scalar field may indeed increase light 
deflection is to consider it as a dark matter fluid, without any coupling to baryonic matter 
and assuming GR as the correct theory of gravity. If the scalar field is massive enough, it 
will obviously cluster and contribute to the total gravitational potential deflecting light. 

Actually, a violation of inequality fl3.5p is a priori also possible when the scalar's stress- 
energy contributes significantly to the curvature of the Einstein metric g^^. Indeed, specific 
models could predict a negligible contribution to the time component g^Q but a significant 
one to the spatial components g*y In such a case, test particles would feel the standard 
Newtonian potential — 1(1 + goo)c^ ~ —GM/r [and the extra scalar force caused by the 
conformal factor A^l^p), which may be tuned to be negligible], but light bending would be 
directly affected by the scalar-induced corrections entering g*j. Therefore, Mll^^ > M^™* 
seems possible. However, such models may be difficult to construct if one imposes the 
positivity of energy, as in Ref. 2^. We will not attempt to study them here, since there is 
anyway no ex perimental evi dence for Mlll'^ > Ml°^. Weak-lensing observations only impose 

Mix ~ A^tot [ni, [ni, E23 • 

Let us end this section by recalling that the conclusion of Ref. 29j] remains correct in a 



limiting case which is often implicitly assumed in MOND-like field theories [42, Il2l[ |: If the 
stress-energy of the scalar field is negligible with respect to that of baryonic matter, then 
the curvature of the Einstein metric (7*^ is basically generated by Mb alone. Since light only 
feels (7*^, weak lensing provides a measure of ~ In other words, light is only 

deflected by baryonic matter, whereas experiment tells us that the much greater amount 
of "dark matter" does deflect it too. We will review in Sees. IIII CI and IIII Dl how the same 
authors as Ref. 2^ devised a clever way to solve this crucial difficulty of RAQUAL models. 



B. A RAQUAL example increasing light deflection 

Let us illustrate that RAQUAL models may predict a much larger light deflection than 
if there were baryonic matter alone in GR, i.e., that ^ Mb is possible. We consider 

the action 

S = I ^V^.{r* - C-'h-] + scatter ~g,u = A'{^)g;^] , (3.6) 

where as before s = g'^^d^ipdyip, io is a length scale (e.g. io oc c^/ao in MOND-like models), 
and n is a positive odd integer. As underlined at the end of Sec. [Tll such a model satisfies 
the two consistency conditions (a) and (b), except on hypersurfaces where s vanishes. It 
suffices to add a small standard kinetic term —es (with e > 0) to the above action to cure the 
singularity on such hypersurfaces. However, we will not take into account this refinement 
because the above model will anyway not reproduce the correct MOND dynamics. 
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The field equations deriving from action (13. 6p read 



- = ^ (^,:r"" + T-^) . (3.7a) 

„4"-'v; (."-'ar*,) = (3,7b) 



wliere tlie stress-energy tensor of tfie scalar field is given by 



^4 



= (^^ s)''-' [nd.^d,^ - \sgiy^ . (3.8) 



Let us clioose for instance an exponential (Brans-Dicke-like) matter-scalar coupling function 
A{ip) = e""^, so that d\n A{ip) / dip = a is constant. The solution of Eq. (]3.7bp . around a 
spherical body of mass M, is then easy to write at lowest order (i.e., neglecting the curvature 
of spacetime): 

^ = ^°+2;73^l^^J V ^'-'^ 

The stress-energy tensor of the scalar field reads thus 



.4 £ 2^2 y ^^2^2 / \ M ^^M^ 



and Einstein's equations (I3.7ap can now be solved straightforwardly, still to lowest order 
around a flat Minkowski background. In particular, one may write {87iG/c'^)(Tqq —^q^qT*) = 
Rqq ~ iAln(— (jfQQ), and this gives immediately 



rc2 n(2n-3) V ^^^c^ J VcV ' 

The i?*^ component of Einstein's equations can now be used to solve for g*^., but one must 
specify the coordinate system to write the solution explicitly. In order to compute easily 
the light deflection angle below, it will be convenient to use the so-called "Schwarzschild 
coordinates" (which define an "area radius" r), such that the line element takes the form 

dsl = -B{ryde + A{r)dr^ + r\de^ + sin^ 9 d^^). (3.12) 

One then finds straightforwardly 

The physical metric g^u = e^"'^^'*^, to which matter is coupled can now be deduced from 
Eqs. (13.91) . (13. lip and (I3.13p . In particular, the gravitational potential felt by test masses 
reads V = — 1(1 + goo)c^, and its radial derivative drV gives their acceleration 



GM f2aGMio^ i/(2n-i) 



+ 



ac^ f 1\ f 2 \ aGM 



(3.14) 
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For large enough values of the integer n, this model predicts thus a MOND-like acceleration 
oc 1/r, together with the standard Newtonian one GMjr^ and a constant contribution 
ac^/d-Q. This constant contribution comes from the conformal factor A^{lp) entering the 
physical metric, i.e., from a scalar exchange between the massive source M and the test 
mass, as in line (a) of Fig. O above. The standard Newtonian force comes from the curvature 
of the Einstein metric g*^^ caused by the matter stress-energy tensor. On the other hand, 
the 1/r force comes from the curvature of (yf*^ caused by the stress-energy tensor of the 
scalar field, as illustrated in line (c) of Fig. [Si However, let us underline that this model 
does not reproduce the MOND dynamics, first because the constant contribution happens 
to dominate over the 1/r force, but above all because the coefficient of this 1/r force is 
proportional to the baryonic mass M instead of a/M. 

In order to compute t he lig ht deflection angle predicted by such a model, let us use its 



exact integral expression [122 



-7r + 2 



dr 



A{r)B{r) 



1/2 



(3.15) 



valid for any metric expressed in Schwarzschild coordinates f l3.12p . In this integral fl3.15|) . 
ro denotes the smallest distance between the light ray and the center of the deflecting body, 
which may also be replaced by the impact parameter h = tq + 0{l/(?) at lowest order. 
Actually, since the coordinates we chose to write Eq. fl3.13|) are such that the Einstein 
metric g*^^, takes the Schwarzschild form, this is no longer the case for the physical metric 
g^y = A'^{ip)g*^j^, because the radial dependence of the scalar field (13. 9p spoils the exact 
expressions gg^ = r^ and g"^^ = r^sin^6'. However, we saw in Sec. IIII Al that light is totally 
insensitive to any global factor of the metric, so that one may actually compute its deflection 
directly from the Einstein metric (7*^. The anomalous terms in Eqs. (I3.11|) and (I3.13P are of 
the form 

^(r) = l + ^, S(r) = l + ^, (3.16) 

where p = 2/(2n — 1) = 0{l/n) is a small number, and where = 0{l/n) whereas 
= 0{n). For any value of these coefficients, integral fl3.15|) gives straightforwardly 



p Tip/2) 



and in our present case 



bP 



AGM TiaGM /2aGMio 



6c2 



l/(2n-l) r 



1 + ( - 

n 



4 



(3.17) 

(3.18a) 
(3.18b) 



It is interesting to note that in Eq. (I3.18ap . pk^ = 0{'nP) dominates over /c_4 = 0{n^^), in 
spite of its small factor p. For large values of n, the anomalous contribution iraGM / Iqc^ 
to the deflection angle (]3.18bp is precisely half the result GR would have predicted in 
presence of pressureless dark matter reproducing the 1/r acceleration in (13.140 [namely, 
Tqq^^ = aMc^ /Anior'^ and Tj^^^^ = 0]. This is due to the specific form of the scalar stress- 
energy tensor (I3.10p . whose spatial components T*^ do not vanish. But besides this subtle 
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factor 2 discrepancy with the standard dark matter prediction, let us also recall that light is 
insensitive to the conformal factor responsible for the dominant anomalous contribu- 

tion in fl3.14p . namely the constant one ac^/£o- In conclusion, using the notation introduced 
in Sec. IIII At the present model actually predicts M^^"^ -C -Mt™*, which is even worse than 
the inequality fl3.5p derived in standard scalar-tensor theories. On the other hand, for a large 
enough impact parameter b, the light deflection angle (]3.18bp . A6 ^ naGM/ioc^ = const, is 



much larger than the one predicted by GR in presence of M alone, namely AGM/bc^ — > 0. 



Therefore, we have exhibited a RAQUAL model predicting MlH'^ ^ Mb. 

One might try to refine the model (13.61) by choosing another matter-scalar coupling func- 
tion, for instance 

A{^) = exp(a(^-2"), (3.19) 

with the same integer n as in the kinetic term of the scalar field. A first useful consequence 
of such a choice is that the extra constant acceleration we found in Eq. (I3.14p above becomes 
now negligible. Indeed, it came from the conformal factor A'^{ip) relating the physical and 
Einstein metrics, so that this extra acceleration derives now from a(y9~^"c^ = (9(r~^"'), and 
tends quickly towards as r grows. On the other hand, the 1/r acceleration of Eq. (I3.14p 
is still predicted, since it comes from the curvature of the Einstein metric caused by the 
scalar's stress-energy tensor oc For large enough impact parameters, such a model 

predicts thus 

MlT ^ ^M™* » Mb, (3.20) 

the factor ^ coming from the same calculation as Eq. (13.180 above. The second conse- 
quence of the choice (I3.19P is that Eq. fl3.7bp now involves the non-constant coefficient 
d\nA{ip)/d(p = —2na/(p'^"'~^^. If one naively counts the powers of dr^p and cp entering this 
Eq. fl3.7bp . one could then deduce that y^^" oc T^atter °^ Using now Eq. (13.81) . this would 
imply T*"^ oc {dr^Y^ oc v^M, which looks like the dark matter MOND wishes to mimic. 
However, this would be the case only if 99 — >■ as r — > 0, within body M. Then one can 
prove that if oc ^/r for small radii, and J d?'xT*'^ would diverge within the body. Therefore, 
such a solution is unphysical. On the other hand, if y^o 7^ within the body, then one 
gets back an exterior solution of the form (13. 9p . where a is replaced here by —2na/Lp'^^^^ so 
that the anomalous force is a priori again proportional to M instead of the looked-for \fM. 
There could still remain a possibility to recover the MOND force with such a RAQUAL 
model, if one could prove that <^ — ^ <^o oc M^/^*^ within the body. This would be the case 
for instance if ^ at a given radius, say at the surface of the body. However, this 
would correspond to a diverging matter-scalar coupling strength at this radius (singular and 
obviously forbidden by solar-system tests), and we anyway do not see how to impose such 
a relation, nor the precise value of ipQ. In standard scalar-tensor theories, the scalar field 
tends towards a constant ipo when r — > 00, this constant being imposed by the cosmologi- 
cal evolution of the Universe, and being constrained by post-Newtonian tests of gravity in 
the solar system. Here, the scalar field = y^o + const, r^"^/'^^""^-' diverges at infinity, as 
soon as a single particle exists in the Universe. This pathological behavior is linked with 
the fact that the solution of Eq. (13.7bp cannot be written unambiguously, the integration 



constant <^o remaining undetermined. In the model (13.90 . where A{(p) = e°"^, this constant 
did not enter any physical observable, but it does when the matter-scalar coupling function 
is no longer a mere exponential, as in our present attempt (I3.19p . Therefore, although this 
matter-scalar coupling function looks promising to reproduce the MOND dynamics from a 
simple RAQUAL model, it fails to define a fully predictive theory. Its only solid prediction 
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is that light bending would be twice smaller than the prediction of GR in presence of dark 
matter, c f . Eg . (I3.20p . As far as we know, this is still allowed by weak-lensing observations 

(nsl. fnol. [T2of 



In conclusion, we have exhibited RAQUAL models predicting a much larger light deflec- 
tion than GR in presence of baryonic matter alone, but this framework does not seem to 
reproduce the MOND dynamics consistently. Note that the class of RAQUAL theories ex- 
amined in the present subsection IIII Bl are actually models of dark matter. Extra Newtonian 
force and light deflection are caused by the stress-energy tensor of the scalar field, which 
plays thus the role of dark matter. The only difference with standard dark matter models 
is that its distribution is strictly imposed by the one of baryonic matter, whereas they are 
usually only related via the dynamical evolution of large-scale structures. 



C. Disformal coupling 

A nice trick exists to increase light deflection within scalar -ten sor theories or RAQUAL 
models. It is similar to Ni's "stratified theory of gravity" [H, Il23 |. and was introduced by 



Bekenstein in [28|; see also |29|. Instead of assuming that matter is coupled to a physical 
metric = A'^{ip)g*^j^, conformally related to the Einstein metric g*^^^, one may consider a 
"disformal" relation of the form 

g,,^ = A\<^, s)gl^ + B{ip, s)d^yDdu^, (3.21) 

where as before s = g'^^d^ipd^ip. The crucial difference is that the d^ipdyip contribution now 
privileges a particular direction, namely the radial one in spherically symmetrical systems. 
This extra contribution may remain negligible in the equations of motion of test particles, 
but now light can be directly sensitive to the presence of the scalar field, contrary to the 
case of conformal coupling that we discussed in Sec. IIII A[ 

Let us emphasize that because of the disformal coupling above, light (or matter) can 
travel faster or slower than gravitons depending on the sign of 5. In a similar way as in 
Sec. Ill Et this does not imply a violation of causality, provided that the matter metric is still 
of Lorentzian signature. The matter metric Eq. fl3.2ip is causal if 

A^{^,s)>Q and A^{^, s) + sB{^, s) > Q (3.22) 

for all (p and s. Notice that if A does not depend on s, B must depend on s for the latter 
inequality to hold. Beware that the se k ii id o f models used in the literature as realizations 
of varying speed of light scenarios 124, 125l ] do not respect this last inequality. Indeed, 



whereas A was taken to be 1, B was chosen constant B = —L^ with L being some fixed 
length scale. Therefore matter field equations are elliptic at local scales (s > l/L"^), and do 
not have a well-posed Cauchy problem. 

By analogy with the Schwarzschild metric, whose positive contribution in the radial 
component grr = 1 + 2GM/rc^ + 0{l/c^) contributes positively to light deflection (together 
with the standard Newtonian potential involved in (^oo), one understands intuitively that 
positive values of B{s,(p) should increase light deflection. To prove it, let us first recall 
again that light is insensitive to any global factor of the metric, because of the conformal 
invariance of electromagnetism in 4 dimensions, Eq. fl3.ip . Therefore, instead of the metric 
fl3.2ip . one may equivalently consider g^i, = g^^ + {B / A?)dfj^(pdy(p. Let us then assume that 
the coordinates have been chosen so that the Einstein metric g*^^^ takes the Schwarzschild 
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form (13.121) . Since, in a static and spherically symmetric situation, the disformal contribution 
{B / A^)d^ipdyip only contributes to the radial component (]„, the metric Qfj^y is also of the 
Schwarzschild form, and one may thus use integral (I3.15P to compute the light deflection 
angle. The only difference with GR is that the radial metric component now reads A = 
>^GR + {B/A'^){drfy. [Beware not to confuse the metric components A = (jrr and B = —goo 
defined in (13.121) with the scalar-field functionals A{s, (p) and B{s, (p) entering (I3.2ip .] At first 
post-Newtonian order, all contributions coming from g*^^ reproduce the general relativistic 
result, and the only extra contribution is thus the one proportional to {B/A'^){dr'p)'^, assumed 
to be of order 0{l/c^) too: 



Even without performing the explicit integration in particular models, it is thus clear that 
a positive value of B{s,ip) always gives a positive contribution to light bending, since the 
extra term in (13.231) is the integral of a positive function. 

For instance, if one tunes the model such that {B/A'^){drfy = 4:\/GMao/c^ is indepen- 
dent of r (as we will do in Sec. IVl below), then 



and the second term reproduces exactly the deflection GR would have predicted in pres- 
ence of dark matter (which is independent of the impact parameter, as if there were a 
deficit angle in a flat conical space). Indeed, it is straightforward to prove that a dark 
matter density Tqo/c^ = V GMao / ATrGr"^ generates a logarithmic MOND-like potential in 
the time component of the metric, — ^oo = 1 — 2GM/rc^ + 2\/GMao Inr/c^ + O^l/c^), and 
that the corresponding radial metric reads (]„■ = 1 + 2GM/rc^ + I^JGMaoj (? + C(l/c^) 
in Schwarzschild coordinates. A similar calculation as above then yields the same light 
deflection angle (I3.24p [see Eq. (I3.32p below for an explicit derivation]. Note in passing 
that g„ 7^ — l/^oo iii such coordinates, contrary to the standard Schwarzschild solution 
in vacuum. On the other hand, if one rewrites this same metric in isotropic coordinates, 
(is^ = —gooc^dt"^ + (jpp [dp^ + p^dO"^ + sin^ 6'(i</)^) , one indeed recovers the usual relation 



gpp = -l/goo + 0{l/c^) = 1 + 2GM/pc^ - 2^GMao In p/c^ + C(l/c^). This remark reminds 



us that a metric is coordinate-dependent, and that observable quantities cannot be read 
naively in its coefficients. Even the signs of the predicted effects are not obvious. For in- 
stance, if one had naively tried to impose A'^{s, ip) = {l + 2^/GMaohlr /c^) but grr = — l/^oo 
in Eq. (13.210 . then one would have needed a negative B{s,(p), and the predicted deflection 
angle would have been negative (and diverging), A9 = A6'gr — 2n\/ GMao In 2h/(?. 

Equation (I3.23P proves that a positive value of B{s, ip) always increases light deflection. 
However, hght rays satisfy then gp^^dx^dx"^ = g^^dx^^dx"^ = —B/A'^{dp(pdx'^y < 0, and 
photons are thus timelike with respect to the Einstein metric g*. This has been considered 



Note that although the deflection angle is constant, the direction of deflection is radial with respect to the 
deflecting body, so that distant objects are deformed and weak lensing observations can be performed. A 
constant deflection angle in a constant direction would have just shifted globally the image of the distant 
objects, without any deformation. 




(3.23) 




(3.24) 
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as a serious flaw by the authors of Refs. 
perturbations of the Einstein metric g'^iA. 
imphed causahty violations. Reference [2 



M, Indeed, this means that gravitons (i.e., 
are superluminal, and it was believed that this 



even presents the reasoning the other way: If 
one imposes from the beginning that no field (even gravity) can propagate faster than light, 
then a disformal coupling 03.211) needs -B < 0, therefore light deflection is even smaller that 
the one predicted in standard scalar-tensor or RAQUAL models, see Sec. IIII A[ This sufficed 
for the authors of these references to discard such models. Actually, as already discussed in 
Sec. Ill Et such a superluminal propagation does not threaten causality. This is even clearer 
in the present case, since one is used to write the field equations in terms of the Einstein 
metric g^^^ (i.e., in what is usually called the "Einstein frame"), in order to diagonalize the 
kinetic terms and prove that the Cauchy problem is well posed. To guarantee causality, like 
in GR, one needs to assume that spacetime does not involve any closed timelike curve with 
respect to g^^. Since the causal cone of matter, defined by the physical metric g^j,^, is always 
interior to the one defined by g^^ if i? > 0, there cannot exist any CTC for any matter field. 
In conclusion, the nice trick of disformal coupling, Eq. 03.211) . is consistent from the point 
of view of causality. 

It remains however to study its consistency within matter. Indeed, the physical metric 
g^i, is involved in many places of the matter action, to contract indices, to define covariant 
derivatives, and in the volume element \/—g df'x. One should thus take into account the 
matter contribution 5S'matter['V^; ^/xi^]/^'/' in the Euler-Lagrange equation 02.141) for and 
study simultaneously the matter field equations 5Sras.tter[i^',giiu\/5''P = (where as before V' 
denotes any matter field). One should then identify any second derivative of (p which appears 
in them, and prove that the full set of field equations, for matter, the scalar field and the 
Einstein metric g^^, has a well-posed Cauchy problem. This is usually done by diagonalizing 
the kinetic terms, i.e., by separating all second derivatives in the left-hand side of the field 
equations, and by proving that their differential operator is hyperbolic. In the general case 
of a theory involving any kind of matter, this is quite involved and we postpone this crucial 
study to a future work. Indeed, if one considers an action of the form 02.131) with a physical 
metric 03.2ip . the scalar field equation reads 
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(3.25) 



where T^^ = (2c/ \/—g) S Scatter /Sg^u is the physical stress-energy tensor, T = T^^g^y, and 
~g = A\l + sB/A^)g^ from definition ([3211). Note that Eq. ([325]) reduces to Eq. (l2T5il for 
A = e°"^ and 5 = 0. The right hand side of the above equation does not involve any second 
derivative and is merely a source term [if OA/ dp were vanishing, d^p = would always be 
a solution]. The left hand side however contains second derivatives of the scalar field and 
of matter fields, through terms like dT. The equation is not diagonal and, moreover, the 
effective metric multiplying the second derivative of p has a rather complicated expression 
that notably involves terms like T^'^, T^^^Oppd^p. The Lorentzian signature of the latter is 
therefore a quite involved question, even if we expect that it may depend on generic energy 
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conditions in the matter sector. Note that the hterature on such MOND-hke field theories 
did not discuss the consistency of the field equations within matter. We will not examine 
either the fully general case in the present paper, but we wish to underline that this should 
be done. We will however analyze below the particular case where matter is described by a 
pressureless and isentropic perfect fiuid. Although much simpler than the general case, it is 
relevant for most matter in the universe, and exhibits the subtle conditions which must be 
imposed on the coupling constants A and B for the scalar field equation to remain always 
consistent. As far as we are aware, this is the first time such a discussion is presented. 
The action defining a pressureless and isentropic perfect fiuid reads ^matter = 

— J d^x^—g pc, where p is the conserved matter density as observed in the Jordan frame. 
If = dx^ / -sj —g^j,i,dx^^dx^ denotes the matter unit 4-velocity, one has thus V\{pu^) = 0, 
which can be equivalently written in terms of a partial derivative as ^\{^J—g pu^) = 0. 
This equation implies that the matter current density ^/—g pu^ is determined everywhere 
once initial conditions are given at a point of the fiowlines, and it is therefore unchanged 
if the metric g^y is varied: pu^) / 5g^u = 0. In particular, the coordinate-conserved 
density p = \/^-g pni^ is independent of the Jordan metric g^y, and thereby indepen- 
dent of both the Einstein metric (7*^ and the scalar field (y9 (as well as its derivatives dy^^p). 
The pressureless perfect fiuid action may thus be written as S'matter = — f d^xpc/vP = 

— / d'^x pc^y—g^iydx'^dx'', which does depend explicitly on the scalar field via the Jordan 
metric fl3.2ip . and therefore contributes to the scalar-field dynamics within matter. Intro- 
ducing the Einstein-frame unit 4-velocity = dx^ / ■^—gj^j^dxf^dx" (independent of and 
its derivatives), one may thus write the action of the scalar field as 

(3.26) 

where all the scalar-field dependences are now explicitly exhibited (neither y/—g^ nor pc^/u^ 
depend on it). In conclusion, the dynamics of the scalar field within matter is described by 
an action generalizing the RAQUAL form fl2.13p . in which the function f{ip, s) is replaced 
by 

/(^, s, <9^^) = /(^, s) + ^ AJl-{u^.d,ipfBlA\ (3.27) 

where A{ip, s) and B{ip, s) are themselves functions of the scalar field and s = g'i'^d^ipdyip. 
Since the derivatives of the scalar field do not enter only via s but also via u'^d^ip, the 
conditions that such a function / must satisfy to define a consistent scalar-field equation 
(hyperbolicity and Hamiltonian bounded by below) are significantly more complicated than 
(a) and (b) of Sec. Ill El above [that / itself must still satisfy to define a consistent scalar-field 
equation in vacuum]. Using as before a prime to denote partial derivation with respect to s, 
i.e., /' = df{ip, s, u^d^(p)/ds, and introducing the shorthand notation p = 2nGp/ \^—g* m°c^, 
C = Bj ^ A? — {utd^pYB, and D = C + {u'^dfj_(pyC^ / B, the hyperbohcity conditions for 
the scalar field equation within matter may be written as 

(al) /' > 0, 



The coordinate-conserved density p is often denoted as p* in the general relativistic hterature, but this 
star index would be confusing in our present context where Einstein-frame quantities already bear such 
an index. 
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(bl) sf" + f + pD-4p{u':d^if)^C' 

± I (sr - pD) - 4p (u^d.^y (r + 2pC") {2sC' + D)\ > 0. 

Note that in condition (bl), the inequahty must be satisfied for both signs of the square 
root, and it suffices thus to check its lowest value, involving a minus sign. At the lowest 
post-Newtonian level, in the limit (u^d^ip)"^ ^ k|, one may simplify these conditions as 

(a2) f > 0, 

(b2) 2s f" + r > 0, 

(c2) f' + 2pC>0. 

If B is chosen positive in Eq. fl3.2ip . in order to increase light defiection in the MOND 
regime, then pC > and condition (c2) is implied by (a2). On the other hand, in the 
ultrarelativistic regime, or in the cosmological one where \do(p\ ^ \diip\, conditions (al) and 
(bl) reduce to 

(a3) f > 0, 

(b3) 2s f" + f + 2p {AsC + D)>0. 

Note however that conditions (al) and (bl) should in principle be satisfied in any situation 
(within matter), in order for the model to always remain a well-defined field theory. 

In addition to (al) and (bl), the contribution of the scalar field to the Hamiltonian should 
also be bounded by below. In the rest frame of the perfect fiuid, such that = (1, 0, 0, 0), 
and simultaneously with coordinates which diagonalize locally the Einstein metric (7*^, = 

diag(— 1, 1, 1, 1), this contribution is proportional to / + 2{doip)'^{f' + 2pC). The general 
conditions needed to ensure that it is bounded by below are more involved. It suffices to 
check that this is indeed the case for each particular model that one may consider. 

Let us finally mention that no extra condition is necessary to ensure that the argument 
of the square root is positive in Eqs. fl3.26l) - fl3.27l) . Indeed, since matter is assumed to be 
metrically coupled to the Jordan metric ^^j^, we know that g^^u'^u'^ = —1, i.e., that the 
matter 4- velocity must be timelike with respect to ^^j^. This implies 1 — {u'^dfj^ipY'B / > 0, 
which should be understood as a condition (automatically) satisfied by the Einstein-frame 
4- velocity m^, but not as a constraint on the functions A{ip, s) and B{{p, s). 

Besides this important question of the consistency of the scalar field equation within 
matter, it also remains to study the precise prediction s of this class of models. Specific 
examples have been analyzed in Refs. @, H, H, Il25[ |. We will examine a very different 



one in Sec. IV Bl below, but still within this class of disformally coupled scalar-tensor theories. 
We will notably show that the predicted light defiection can be consistent both with weak- 
lensing determinations of dark matter and with post-Newtonian tests in the solar system. 

D. Stratified theories 

Because the above superluminal propagation of gravity was considered as deadly, Sanders 
proposed in 32, 3^ another kind of disformal coupling, where the physical metric to which 
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matter is universally coupled takes the form 



g,, = A\^)g;, + Bi^)U,U,. (3.28) 

This is thus strictly the same form as Eq. fl3.2ip above, but the gradient d^(p has been 
replaced by a vector field (not to be confused with the matter uni tary velocities u'^ and 



^ary 

introduced in Sec. IIII CI above). Bekenstein's theory TeVeS [30|, l3l| also assumes the 
above disformal relation. The dynamics of the vector field is defined by a specific kinetic 
term of the Proca form (9[^f/j^])^, and its norm is imposed to be g'i^UfJJy = —1 thanks to a 
Lagrange parameter. The reason why it is necessary to fix its norm is because the magnitude 
of scalar-field effects would otherwise depend on the background norm of the vector field. 
Since backgrounds where = would be possible, one would not predict any extra light 
deflection caused by the scalar field in such situations. Moreover the Lagrange multiplier 
ensures that the vector field equation is consistent with any type of matter one considers. 
Indeed, without the Lagrange multiplier, the equation of motion reads 



v;(vi^f/:i) oc y-Af'^'^u,, (3.29) 

where = g^'^U^^, and its divergence imposes that the stress-energy tensor of matter must 
satisfy the non-trivial property ^*y{TP^U —g / s/ —g^) = 0, that may be satisfied only by 
very particular type of matter. On the other hand, a term like \{UpU^ + 1) in the action, 
where A is a Lagrange multiplier, gives rise to a term proportional to in the left hand 
side of Eq. (13.291) . whose divergence equals the a priori non vanishing divergence of the right 
hand side, thus giving its value to A. Note finally that the consistency of the field equations 
requires that the above matter metric fl3.28p has a Lorentzian signature. This condition 
reads > and > B for all ip, where we used U^Ul^ = —1. 

One could think that the vector field is chosen timelike (negative norm), contrary to the 
gradient d^cp in the vicinity of clustered matter, in order to preserve isotropy of space which 
is very well tested experimentally. However, this is not the case, since a timelike vector does 
define a preferred frame (i.e., an "ether"), namely the one where it takes the canonical form 
[/^ = (±1,0,0,0). We will briefly discuss the related experimental difficulties in Sec. IIVBI 
below. 

The bonus added by this vector field [/^, with respect to the purely scalar case (13.211) . 
is threefold: (i) it allows us to increase light deflection without necessitating superluminal 
propagation of gravity; (ii) it can also be used to tune a specific kinetic term for the scalar 
field, in order to avoid its own superluminal propagation identified in Sec. Ill Fl above: (iii) the 
dynamics of the various fields is easier to analyze even within matter, and there does not 
seem to exist any generic difficulty with the Cauchy problem. Of course, the third point 
is actually the only useful one, since we know that superluminal propagation does not 
threaten causality provided there always exists a nonva nish i ng e xterior to the union of all 
causal cones, as discussed in Sec. Ill El above and in Refs. |l06 . IOTI ]. However, the real reason 



why this vector was considered in Refs. |30|, |31|, 132|, |33[ is indeed to avoid any superluminal 



propagation. Point (ii) is easy to understand: If the kinetic term of the scalar field is defined 
as {gt^'^ — U^U^)d^ipd^(p, instead of the standard s = g^'^d^^Lpd^ip, then a violation of condition 
(c) of Sec. Ill El only implies that the scalar causal cone is exterior to the one defined by the 
inverse of gi^'^ — Uj^U^^, namely g*^^ + \UpUi, = diag(— |, 1, 1, 1) in the preferred frame and 
in locally inertial coordinates. It becomes thus possible to reproduce the MOND dynamics 
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while avoiding a propagation of (y9 outside the causal cone defined by y*^ = diag(— 1, 1, 1, 1). 
Since point (i) is subtler, we discuss it now in more detail. 

It is actually quite simple to understand why a coupling of matter (including electro- 
magnetism) to the disformal metric fl3.28p suffices to increase light deflection. The basic 
idea is similar to the case of the previous disformal coupling (13.211) : Photons are now di- 
rectly coupled to the scalar fleld, whereas they do not feel any global factor of the metric. 
More precisely, in the preferred frame where t/^ = (±1,0,0,0), and in locally inertial co- 
ordinates such that = diag(— 1, 1, 1, 1) at a given spacetime point, one notices that 
—Ufj^Uu = diag(— 1, 0, 0, 0) behaves like the ygo component of the metric, and conversely 
that y*^ + U^Ui, = diag(0, 1, 1, 1) behaves like its spatial components g*y The vector fleld 
is thus a mere tool to separate by hand the space and time components of a tensor (im- 
plying obviously a violation of Lorentz invariance, namely the existence of an ether). Its 
interest is that one may now multiply g^^ and g*j by different functions of the scalar fleld. 



Inverse factors e"^""^ and e~'^°"^ have been chosen in Refs. [30|, l3l|, l32|, |33|, in order to mimic 
GR in isotropic coordinates. One could consider more generally inverse factors A'^{if) and 
((/)), but an exponential (Brans-Dicke-like) coupling A{ip) = e""^ suffices to predict the 
right phenomenology. The physical metric is thus chosen as 

yoo = e2"^ySo; = e-'"^^?;,. (3.30) 
This can be rewritten covariantly thanks to the vector field f/^ defining the preferred frame: 

~g,u = -e2"'^f/^f/, + e-2°'^(y;, + f/^f/,) 

= e-'"^y;, - 2 U^U, smh{2a^). (3.31) 

This matter metric has a Lorentzian signature, as can be easily checked. Beware that our 



notation differs slightly from that of Refs. [30, |31[. In these references, the matter-scalar 
coupling constant a reads a/ k/Air where k is a dimensionless parameter. 

In order to reproduce the MOND dynamics at large distances, the model must be tuned 
so that e^"*^ = 1 + 2 y/GMoo In r/c^ in this regime. Moreover, if one insists on the fact that 
photons should propagate faster than gravitons (even if it is actually not necessary), one 
must have aip > at every spacetime point. The logarithm entering this MOND potential 
should therefore be understood as ln(r/£), where £ is a length scale significantly smaller 
than the size of the considered galaxy (either imposed dynamically by the model, or simply 
a universal constant much smaller than the size of any galaxy). Indeed, ds"^ = dsl = 
2e'^°"^ smh.{2aip){U^dx'^)'^ > because aip > 0, and the cone defined by g^^ is thus wider 
than the one defined by y*^. 

Note that the character spacelike of d^ip or timelike of does not change the sign 
of the square {d^^pdx^Y {U^dx^Y. The crucial difference is that we now consider a 
negative factor for the term U^Up entering the disformal metric fl3.3ip . whereas we needed 
a positive value of B{s,ip) to increase light bending with metric (13.211) . In fact, it may 
seem paradoxical to increase light bending with the negative term entering Eq. (13.311) . since 
photons now propagate faster than in GR (i.e., faster than gravitons), if aip > 0. Intuitively, 
the trajectories of fast particles are almost straight lines, whereas slower ones are significantly 
curved by the Newtonian potential. This reasoning just happens to be erroneous, because 
light bending depends very differently on the temporal (yoo) and spatial {(jij) components of 
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As an illustration of this point, let us recall that the scalar field surrounding a body of mass M is not 
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the metric, as illustrated by Eq. (13.150 in Schwarzschild coordinates. At linear order, Ref. 2^ 
actually proves that this integral involves the radial derivative drQoo of the time component, 
thanks to a partial integration, whereas the radial component c/rr enters it without any 
derivation. In other words, one needs an additional positive term in g^r to increase light 
bending, as proven in Eq. (13.231) above, but one needs a term whose radial derivative is 
negative if it enters ^oo; in Eq. (I3.3ip above. This is the reason why the Newtonian and 
MOND potentials enter so differently in the metric predicted by GR in presence of dark 
matter ^oo = -1 + 2GM/rc^ - 2y/GMao \nr/c^ and grr = l + 2GM/rc^ + 2yJGMao/c^ [see 
below Eq. (13.241) ]. In spite of their relatives signs, all the above terms contribute positively 
to light bending, because those entering g„ are positive, and because the derivatives of those 
entering ^qo are negative. 

The above discussion illustrates once again that observable quantities, or even their signs, 
can be subtle to identify in the components of a metric tensor. The only way to be sure of an 
observable prediction is to compute it explicitly. We already mentioned, below Eq. (13.241) . 
that the GR metric in presence of dark matter takes precisely the form (13.301) in isotropic 
coordinates. Let us however derive now the light deflection angle it predicts. First of all, 
since light is insensitive to any global factor, let us multiply it by e^°"^, i.e., consider instead 
the metric g^Q = e^^^^g^Q and (jij = g*j. If the coordinates have been chosen so that the 
Einstein metric g'^^, takes the Schwarzschild form (13.121) . this is also the case for ^^i, and one 
may thus use the integral expression (I3.15p . At first post-Newtonian order, all contributions 
coming from (7*^ reproduce the general relativistic result, and the only extra contribution 
is proportional to the scalar field, namely aip = aipo 
MOND regime. One thus gets straightforwardly 
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(3.32) 



which confirms that the light deflection angle is positive although photons are faster than 
in GR (with respect to the Einstein metric g^i,)- This angle is exactly the one predicted 



only given by (^locai = —aGM/rc^ (or \/ GMuq hir/ac^ , depending on the regime), but more generally 
by V'o + <y5iocai when one takes into account the asymptotic value of the scalar field ipo which evolves with 
cosmic time. Since the "speed" of photons (or, more precisely, their character timelike, null or spacelike 
with respect to the Einstein metric g*^) depends on the sign of ipo + ipiocai, one may expect that light 

bending would be affected by (po. Actually this is not the case, as shown by the first line of Eq. (|3.32p . 
This also explains why the time component fcg enters multiplied by a factor — p in Eq. (|3.18ap . contrary 

to the radial component fc^. 

Actually, there is a subtle modification of the Einstein metric 17*^ itself in TeVeS, because of the contri- 
bution of the energy-momentum tensor of the unit vector Up,; see Eq. (58) of |30| and Eq. (6.2) of 31|. 
However, this correction is chosen to be negligible in TeVeS. 
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by GR, when a spherical halo of dark matter is assumed to generate a Newtonian potential 
\/GMaQ Inr. Moreover, the first line makes it clear that the background value of does 
not contribute to the integral. It must be stressed that the sign of ip is not crucial in 
itself; what is important is that its derivative is positive. It is therefore possible to have a 
cosmological background such that a^p < locally whereas the amount of light deflection 
remains the same. This remark will be of great importance in Sec. IIVCI since we will show 
that an empirical fact strongly supports that gravitons should propagate at a speed greater 
(or equal) than the speed of photons (i.e., that aip should be negative). 

Although the above construction seems natural (apart from the existence of a preferred 
frame), let us however underline that it is actually quite ad hoc. Indeed, the relation fl3.30p 
has been imposed by hand precisely to mimic GR's metric, but any other relation between 
(7oo and g*j would have predicted another light deflection angle. Even without considering 
the general case f l3.28p . where ^(v^) and B{ip) are fully independent functions of the scalar 
field, it is instructive to analyze a disformal coupling slightly generalizing Eq. (13.311) . namely 

~g,u = -e^^'^U^U, + e-'^^^^ig;, + U^U,), (3.33) 

where k and the matter-scalar coupling constant a are dimensionless numbers to be con- 
strained by experiments. The tightest experimental constraints come from solar-system 
tests, i.e., in the Newtonian regime of such theories, where one must have (f = ipQ—aGM/r(?. 
A nonvanishing background value ipQ of the scalar field is necessary to respect the inequal- 

,2^1 
2 



ity aip > 0, imposed in Ref. [30|] to avoid superluminal gravitons. Since GM/Rc < 



for any body of mass M and radius R, it suffices to impose ipo/a > |. (Let us however 
recall that superlumi nal g ravit ons are actually not a problem, as discussed in Sec. Ill El 
above and in Refs. ^M, [1Q3.) In this Newtonian regime, metric (13.331) reads thus 



^00 = exp(2a(^o) [-1 + 2(1 + a^)GM/rc^] + C(l/c^) and gij = exp{-2Ka(po)5ij[l + 2(1 + 
K,a'^)GM/rc'^] + 0{l/c^), where the constant exponential factors may be eliminated by a 
redefinition of (physical) time and length units. The Eddington post-Newtonian parameter 
^PPN -g ^]-^yg given by 

PPN l + na'^ 1 , {i^-lW /oo.N 
7 = -r- — ^ = 1 + — r- — —■ ('^•'^4) 
1 + 1 + 

The Cassini constraint 901] on this parameter imposes therefore ]k — 1]q;^/ (1+a^) < 2 x 10^^ 



In conclusion, experiment forces us to fine tune n ^ 1, otherwise the matter-scalar coupling 
constant a needs to be small and we get back the serious problem pointed out in Sec. Ill Fl 
above: The MOND logarithmic potential starts manifesting at too small distances, and one 
needs unnatural RAQUAL kinetic terms like the one of Fig. [3] to pass solar-system tests. 
[Note that k = 1 ± 0.1 suffices for the Newton-MOND transition to be much smoother 
that in this Figure, allowing a Newtonian regime up to distances ~ 100 AU instead of 
30 AU, but this remains nevertheless a quite fine-tuned behavior.! On the other hand, if 



K = 1 strictly, as chosen in Refs. |30l . |31| . |32| . 1331 ]. there is no longer any post-Newtonian 
solar-system constraint on the matter-scalar coupling constant a, and one may thus a priori 
choose a = 1 and a RAQUAL kinetic term of the natural form of Fig. [2J [The fact that 
K = 1 implies 7^^^"^ = 1 has been noticed in Ref. [13] ■ However, its author did not choose 
a = 1 (i.e., k = An in the TeVeS notation) but a smaller value a 5 x 10~^ (i.e., k ~ 0.03). 
The reason was that he wished to neglect the energy-momentum tensor of the scalar field 
oc a^, but it is actually negligible even when a = 1 since it is of second post-Newtonian 
order 0{l/c^).] 
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This discussion illustrates that several experimental constraints have already been used 



to construct the model proposed in Refs. [30|, IMI, l32|, |33|] . Moreover, one cannot claim that 



it predicts the right light deflection, but rather that it has been tuned in order to do so. A 
vivid way to illustrate the fine tuning hidden in the choice fl3.30p - fl3.3ip is to exhibit different 
metrics of the form (13.280 . reproducing GR's light deflection in the Newtonian regime of the 
solar system but a different one in the MOND regime. For instance, if one writes this 
physical metric as^^ 

~g,u = -h{v)U^U, + fr{^){g% + (3.35) 

then the choice /o(v3) = e^"^ and fr{'^) = e"^'^/*^"^"^^'^-' reproduces the above model when 
-C but predicts a smaller light deflection when ip > 1/e, notably when oc Inr — oo 
in the MOND regime. On the contrary, the choice fo{^) = e^*^ and fr{^) = e~'^'^^^~^^'^^ 
would also mimic the above model for ^ but predict a much larger light deflection 
when ip > 1/e, in the MOND regime. Therefore the speciflc choice fl3.30p - fl3.3ip . mimicking 
GR, may almost be considered as a flt of some experimental data. It is nevertheless rather 
natural, and the fact that it allows to reproduce GR's predictions in presence of dark matter 
is anyway a great achievement, since light deflection was a crucial difficulty of previous 
constructions. 

Further analysis of such stratified theories however undermines this conclusion, and we 
will see that they cannot survive other phenomenological constraints, unless they are much 
more fine-tuned than they already are. 



IV. DIFFICULTIES WITH STRATIFIED THEORIES 
A. Stability issues 

One of the great advantages of the stratified TeVeS metric fl3.28p . over the disformally 
coupled scalar field fl3.2ip . is the simplicity of the field equations within matter. However, 
as we have seen, the vector field needs to have a constant norm unless consistency problems 
occur in its field equation, and this has two serious drawbacks. First, this leads to a mixing 
of the vector and the spin-2 degrees of freedom and thereby to a nondiagonal system of 
differential equations. Th is means that the actual degrees of freedom have to be computed 



carefully (as done in [11 11 ] in the specific case = g^^). Second, the fact that the vector 
field must be of constant norm spoils its stability. Indeed, as was shown in Ref. jl^], the 
Hamiltonian of the vector field is not bounded by below, and the theory is unstable. The 
dynamics of the vector field derives from the action 

^d'x = J ^d'x [g^/g:'^F^,F,„ - 2\{g^^U,U, + 1)] , (4.1) 

where the stress Ff^i, is given by d^Uu — d^U^, C is the Lagrangian density, and is a 
positive dimensionless number. The Hamiltonian can be computed in the standard way, 
noticing that both A and Uq must be treated as Lagrange multipliers since the action does 



To simplify the discussion, we assume here that the asymptotic value (po vanishes, or, stated differently, 
we merely denote here Lp — Lp^ as ip. 
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not involve their time derivatives. The momenta read then 



where latin indices run from 1 to 3, and indices are raised using the inverse Einstein metric 
gl^'^. In the limit of decoupling gravity ((7*^^ —>■ 1]^^), the Hamiltonian reads 



H = j d^x {71'doUi - C) 



where we used that U^U^ = — 1 which derives from the equation 5S/5\ = 0. One then writes 
doUi as Foi + diUo and integrates by parts the last term. The Hamiltonian, as a functional 
of the independent field configurations Ui and vr*, then reads 



H 



[ d'x l^Tz' + ^(V X Vf + v/TTTP Vtt 
J [Kc^ AttG 



(4.4) 



where bold-faced symbol denotes three-vectors. The first two terms are just proportional to 
^2_|_2j2 Tffi^QYQ Qj^(^ Q^^Q usual electric and magnetic fields associated to U^. Contrary 
to the Maxwell case where Vtt oc VE = in vacuum, the last term does not vanish here 
(even in vacuum), and can take arbitrarily large and negative values. This can be made more 
explicit, for instance, by considering the special case where the vector field derives from a 
static potential 0, i.e., Ui = 9j0(a;); see Ref. jl^]. The Hamiltonian is thus unbounded by 
below, and the theory is unstable, at least as it is defined in action (14. ip . Note that empty 
space {Ui = 0) is not even a locally stable vacuum, because when U. = di(f){x) 0, the 
dominant term is the last one in Hamiltonian (14.41) . Of course, there always remains the 
possibility that higher-order terms could stabilize the model (around a different vacuum), 
but one should exhibit such corrections explicitly to prove the consistency of the model. 

We stressed in Sec. Ill Dl that a competing model (phase coupling gravity) also reproduces 
the MOND dynamics. The theory is however unstable because of the presence of a negative 
potential. Recently, an improved version of this theory (called BSTV) was given in Ref. (33| . 
It involves, besides two scalar fields, a dynamical unit vector field in order to reproduce the 
right amount of light deflection using a disformal metric of type fl3.3ip . This improved 
theory is however still unstable. Three different terms of the full Hamiltonian can take it 
to arbitrary large and negative values. To begin with, there is the instability driven by the 
vector field that we described above, if its kinetic term if of the Proca form (like in TeVeS). 
The dynamics of the two scalar fields q and ip derives from the action^*^ 



S=-l V^d'x Q(5m?)' + [Hq)9r + (/(g) - h{q))U!:U:]d,ipd,ip + V{q) 



(4.5) 

Because the author of Ref. |33|] wrote a quadratic potential V{q) = —Aq^ with A > 0, the 
field q is tachyonic and thus unstable. Note that this potential is not crucial to obtain the 



18 



Beware that Ref. |33j uses the opposite (mostly — ) signature. We have translated its action in our mostly 
convention. 
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MOND dynamics, but is merely added to produce oscillations of the scalar field that may 
generate the cosmic microwave background (CMB) spectrum predicted in cold dark matter 
theories. Curing the instability by changing the sign of A may however spoil this interesting 
phenomenology. The computation of the full Hamiltonian reveals another instability. The 
momenta read Pg = doq and = 2hdoip — 2{U^d^ip){h — f)Uo. This last equality can be 
inverted to find doip as a function of (notably) p^. Since the resulting Hamiltonian has a 
quite complicated expression, we only report its value on the particular field configuration 
such that Ui = 0. Up to the contribution of the vector field, we have: 

Changing the sign of A would not be sufficient to guarantee the stability of the theory. 
Indeed, the functions h and / must behave as h{q) = and /(g) = q^ (at least in some 
regime) in order to obtain MOND dynamics. Therefore the term p^/4(2/i — /) can be made 
arbitrarily large and negative. 



B. Preferred-frame effects 



Ni's stratified theory [123 



is a particular case of disformal coupling (I3.3ip where the 
Einstein metric is no longer dynamical, and assumed to be flat: (yf*^ = Although it 
was huilt to reproduce the Schwarzschild solution at linear order, by imposing e^'^'^ = 1 — 
2GM/rc^^ it is nevertheless known to be inconsistent with preferred-frame tests. Indeed, it 
predicts a post-Newtonian parameter ai ~ —8, whereas the present experimental constraint 
is |ai| < lO"'^ [H, [2|. One may thus wonder if the more general models of Sec. IIII Dl present 
the same difficulty. 



Their experimental constraints have been discussed in Ref. [32j, as well in Ref. fl26l 



m 



a more general context. In the case where the preferred frame defined by f/^ is assumed 
to be the CMB rest frame (as is usually done when analyzing the weak-field predictions of 
vector-tensor theories), Ref. j32| has proven that the post-Newtonian parameter ai ~ — 16a^ 
(where a without any index denotes as before the matter-scalar coupling constant, related 



to the parameter r] of Ref. [32j hy r] = Aa ). The other PPN parameters characterizing 
preferred-frame effects vanish identically, like in general relativity (^2 = as = 0). The 
present experimental bound on ai would thus impose < 6 x 10~^, even tighter than 
the light-deflection constraint < 10~^ on conformally-coupled scalar-tensor theories. The 
fine-tuning problem discussed in Sec. IIIFI above would thus be even more serious in the 
present disformal framework. 



However, the vector field f/^ is dynamical, and Refs. |30l . |31| claim that its field equation 
drives it to be parallel to the matter proper time, i.e., to read t/^ = (±1,0,0,0) in the 
matter rest frame. In such a case, there would be no preferred-frame effect due to the 
motion at 370 km.s~^ of the solar-system with respect to the CMB, but there could still 
exist observable effects due to the planets' velocities around the Sun. The Earth's velocity 
at about 30 km.s~^ would a priori induce effects constrained by the tracking of artificial 



satellites (l27|, and although the bound on might be increased by a factor 10, the fine- 
tuning described in Sec. IIIFI would still be necessary. It remains possible that the vector 
field is also driven to be parallel to the Earth's proper time in its vicinity, but there would 
always exist a transition region between the barycentric and geocentric frames in which test 
particle should undergo observable effects due to the variation of U^,. 
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Anyway, the derivation proposed in [30, [M] actually only proves that = (±1,0,0,0) 
is a solution, but not the solution. In other words, it just tells us that a situation where 
matter is at rest with respect to the preferred frame is possible. This obviously does not prove 
that the model is consistent with high-precision tests of local Lorentz invariance of gravity. 



There remains to perform a careful analysis of the specific models proposed in [30|, l31|, l33 
in order to draw definite conclusions about this problem of preferred-frame effects. We will 
not discuss it any longer in the present paper, because stratified theories anyway present 
several other serious difficulties. 



C. Cherenkov radiation and high-energy cosmic rays 

Throughout this paper, we insisted on the fact that superluminal fields do not neces- 
sarily threaten causality, provided they propagate along cones in spacetime, i.e., if their 
equations of motion remain hyperbolic. We even announced in Sec. Ill Fl that experiment 
strongly suggests that no field should propagate slower than light. Indeed highly relativistic 
matter coupled to a subluminal field emits Cherenkov radiation and thereby loses energy. 
In particular, the various gravitational fields one may consider should not propagate slower 
than light, otherwise ultra high energy cosmic rays would emit "gravi- Cherenkov" radiation 



(following the term used in [35|, l36|), and those that we detect should thus have been pro- 
duced near the Earth. However, no source has been identified within a range of the order of 
the kiloparsec, and this yields a tight bound on how much photons (and hence relativistic 
matter particles) are allowed to travel faster than gravity. 

Cherenkov radiation occurs in electrodynamics when charged particles propagate faster 
than light in a medium with refractive index n > 1. The emission of a real photon through 
the vertex of QED is kinematically allowed in that case because electromagnetic waves do 
not follow the usual dispersion law k^k^ = 0, but satisfy the law |k| = n\kQ\ in a refractive 
medium and have thus a larger wave vector than angular frequency if n > 1; see 35, 3^. 

If light propagates faster than gravitational waves, as is the case in the above stratified 
theory with aip > 0, a similar phenomenon occurs: (ultra relativistic) matter particles 
propagating faster than gravity emit gravitational waves. The same reasoning as above 
shows that this tree process is kinematically allowed. We may however stress a subtlety 
here. Since spacetime is endowed with two different metrics of Lorentzian signature, one 
can define two different sets of locally inertial coordinates (for which one of the two metrics 
(7*^ or Qfj^iy, but not both, is written in its canonical form). There exists therefore two 
(equivalent) formulations of conservation laws. Gravitational waves follow the null geodesic 
of the gravitational metric^^ and satisfy thus the equation g'^^kyk,^ = 0, where fc^ denotes 
their wave four-vector. Their dispersion relation in the matter (Jordan) frame reads thus 
g^^'kf.k^ = {U^kf,yB/[A^{B - A^)], and takes the form 

|k| = |A;o|Vl + 2sinh(2a(^) (4.7) 

in coordinates such that g^i, = rj^^ and simultaneously t/f = (±1,0,0,0) locally (which can 
always be imposed by a local Lorentz transformation keeping g^u = rj^u). Using aip <^ 1, we 



Beware that the spectrum of gravitational waves may be affected by the presence of the vector field. 
How ever, when the kinetic term of the vector field is of the Proca type, there is no such modification; see 



Ref. 



128 | 
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get |k| = |fco|(l + 'iotif). In other words, everything behaves as if gravitons were travelhng 
within a refractive medium of index n 1 + 2aip. The emission of real gravitational waves 
by ultra relativistic matter particles is thus allowed once a(p > 0. The authors of Ref. [3^ 
calculated the rate of energy loss by this process and found: 

dE _ Gp\n - If 



dt K^cii 



ight 



for a particle of momentum p. Here we approximate pure numbers by 1, which is correct in 



view of the exact results derived in j35|. [This result is actually obtained for a scalar particle 
of matter; the case of a proton is more involved but leads to similar formulae; see Ref. (36|]. 
This equation is easily integrated for an ultra- relativistic particle E ^ pcught- Assuming 
that the initial value of the momentum pi is much larger than the final one pf and that n is 
approximatively constant along the particle's path, one finds that it can travel only over a 
distance of order 

^ Gp}{n-ir ^^-^^ 

whose value reads 

2 X 10""^"^ 

TT^m (4.10) 

for p/ ~ 10^'' eV. Therefore, cosmic rays detected on Earth with such a high momenta as 
10^° eV cannot have been produced farther than Lmax- Assuming that these cosmic rays 
have at least a galactic origin, we find L-^ax ~ 10^^ m and therefore n — 1 ^ 10~^^. Note of 
course that gravi-Cherenkov radiation can only be emitted by matter particles propagating 
at a speed greater than Cgravitons = cngu/n, where n — 1 ^ 10^^^. Therefore only ultra- 
relativistic particles with a Lorentz-Fitzgerald factor of order 7 > l/y'2(n — 1) ~ 10^ can 
emit gravitational waves in this way. Assuming that the observed high-energy cosmic rays 
with momenta p ~ 10^° eV are protons, we find that their 7 reads 10^^, so that our analysis 
is self-consistent. 

Since the gravitational field produced by the Milky Way is approximatively MONDian 
in the vicinity of the Sun, we have n — 2a(p ~ 2 GMMUky Way oo In (r /£) /c^ , where £ is a 
fixed length scale and where, numerically, GMniiky WayOo / ~ 10"'''. If one insists on the 
fact that a(f should be positive (so that gravity propagates slower than photons), the only 
way inequality n — 1 < 10^^^ could be satisfied would be to fine tune the length scale i to 
the galactocentric distance of the Sun up to nine decimals! In other words, such a theory 
predicts that high energy cosmic rays should not be seen unless they are produced very close 
to the Earth, or if the Sun were very specifically located within our galaxy. This is clearly 
unacceptable, and aip would anyway become negative for slightly smaller galactocentric 
distances. In other words, gravity must propagate faster than light at least in some regions. 

The only reasonable way to avoid the above fine tuning problem is simply to have a(p < 0. 
In such a case, light propagates slower than gravitons and there is no gravi-Cherenkov 
radiation at all.^° Let us emphasize that it is however in great disagreement with the whole 



Note that one may also consider Cherenkov radiation made of wave s of the vector field. Indeed, depending 
on the sign of the parameters, such "Einstein- aether" waves 128j may propagate slower than light. In 
any case, however, it would be possible to require that aip is negative enough so that light (and hence 
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spirit of the TeVeS theory. Indeed the vector field was precisely introduced to produce the 
right amount of light defiection while avoiding any superluminal field. Actually, the above 
discussion shows that subluminal fields coupled to the matter sector are generically ruled out 
by the existence of ultra high-energy cosmic rays. Therefore MOND-like stratified theories 
reproducing the right light defiection must either be based on the disformal relation fl3.2ip 
(where the vector field assumes the form f/^ = d^ip) with 5 > 0, or on the disformal metric 
Eq. fl3.3ip with aip < 0. In conclusion, light should in any case travel slower than gravitons. 
Let us recall that it does not imply a violation of causality, provided all field equations 
remain hyperbolic. 



D. Binary-pulsar constraints 



The fourth and probably most serious difficulty of stratified theories is their consistency 
with binary-pulsar observations. Indeed, even if the choice k = 1, in Eq. fl3.34p . releases 
the constraint on a arising from precise measurement of the Eddington parameter 7^^^ 
in the solar system, binary-pulsar data a priori impose a tight constraint on it too. The 
literature on binary-pulsar tests has studied wide classes of alternative theories of gravity 
[H. [oil. 93I . 94, 129l . 130 . 13 1| but not the disformal coupling (13.280 . A specific analysis should 
thus be performed, but it goes beyond the scope of the present article. Let us just mention 
here that a strong matter-scalar coupling constant a ~ 1 means that the spin-2 {g*^^ — rj^^) 
and spin-0 (y? — (y^o) fields are of the same order of magnitude in the vicinity of massive bodies, 
i.e., in the Newtonian regime of the model. Near a binary system, these fields cannot be static 
because they are driven by the orbital motion. This causes the emission of gravitational and 
scalar waves (and even waves of the vector field f/^, in the present TeVeS model), which 
extract energy from the system and make the orbit shrink. Because of the specific form of 
the kinetic terms (i.e., of the helicities of the degrees of freedom), the energy loss via spin-2 
waves start at the quadrupolar order (9(l/c^), but it is of the much larger order 0{l/c^) 
for the dipolar spin-0 waves [monopolar scalar waves also exist but cause an energy loss of 
order 0{l/c^) instead of 0{l/c) if the two orbiting bodies are at equilibrium]. The orbit 
should thus shrink much more quickly in the present model than in GR, and t his i s known t o 



be inconsistent with several binary-pulsar observations [9l|, 93, 129l . 1301 . 13ll . 132, 133 
Since this reasoning only depends on the local dynamics of the scalar field, its subtle MOND- 
like behavior at larger distances should not change the conclusion. The scalar waves will 
obviously behave in a non-standard way at distances r ^ a/ GM/ao, but they have anyway 
extracted too much energy from the binary system in the near zone. Note also that the 

-f l3.3ip should not change either the above conclusion. 



peculiar matter-scalar coupling f l3.30p 
It tells us that (f does not feel the matter current and pressure T™^**®'' as in standard scalar- 
tensor theories, but the dominant source remains aTg™'^***^'' oc aMc^ and generates a scalar 
field (f — (fo = 0{aGM/rc^). Therefore, we expect pr esent binary-pulsar data to impose 



4 X 10"^, as in Brans-Dicke theory [l2i, [TsqUTsi 



than the 10 ^ bound imposed by solar system 9ol 



Although this is less constraining 
on conformally-coupled scalar fields 
(I2.i2p - fi2.i3p . this binary-pulsar constraint on the matter-scalar coupling constant is thus a 



matter) propagates slowly enough, to ensure that matter particles do not emit such Cherenkov radiation 
at all. 
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priori also valid for the disformally coupled ones fl3.2ip - fl3.28p .^^ 

Note that this constraint is one order of magnitude smaller than the value ~ 2 x 10~^ 
chosen in TeVeS [sij] (where = k/Air in Bekenstein's notation). This reference needed 
to choose > 1.4 x 10~^ in order to predict negligible MONDian effects at distances 
smaller than Saturn's orbit, assuming a natural RAQUAL kinetic term like the one displayed 
in Fig. [2] above. Since this is inconsistent with binary-pulsar tests, we recover again the 
problem discussed in Sec. IIIFI a needs to be small and we need an unnatural kinetic term 
f{s). In particular, the acceleration cannot be purely Newtonian beyond r ~ 150 AU (if 

= 4 X 10~^), and an approximately constant extra contribution ao must exist between this 
distance and a/ GM/oq ~ 7000 AU, where it becomes of the 1/r MOND form (see Fig. [3]). 
Although this is still allowed experimentally, this behavior seems anyway too fine-tuned. 



E. Discussion 



As Sees. IIIIDI and [IV] above underline, the TeVeS model |30|, l3l|, l33| is certainly the 
best field theory reproducing the MOND dynamics in the literature, but it has still some 
drawbacks. One of them is probably unavoidable: It is fine tuned to obtain the wanted 
phenomenology. A second one seems generic but hopefully curable in some specific models: 
It presents some instabilities. Other theoretical problems should be solved, and have already 
received some answers: It involves discontinuities, too many different fields, and the various 
arbitrary functions it contains look quite complicated. But the most serious difficulties seem 
to be experimental, as discussed in the present section. In particular, it very probably does 
not pass binary-pulsar tests in its present form (or needs to be even more fine-tuned in a 
quite unnatural way). We will also mention in Sec. IVII other generic problems that any 
MOND-like field theory should address. Therefore, in spite of the great achievements of 



Refs. [30, |31|, |33|, it remains interesting to look for other theoretical frameworks able to 
reproduce the MOND dynamics in a consistent way. We will examine a rather new one 
in the following section, although it is still in the class of RAQUAL and disformal models. 
Its original feature is that the theory will be particularly simple in vacuum (outside matter 
sources), namely pure general relativity in a first model, and pure Brans-Dicke theory in a 
second one. However, we will see that the consistency of the field equations within matter 
exhibits subtle but crucial difficulties. 



V. NONMINIMAL METRIC COUPLINGS 
A. Changing the matter spin 2 coupling 

In the present section, we adopt a mixing of the modified-inertia and modified-gravity 
viewpoints discussed in Sec. Ill A] above. Like in the modified-inertia framework, the dynam- 
ics of gravity in vacuum will be defined by the Einstein- Hilbert action (12. ip . On the other 



This discussion is reminiscent of a particular tensor-bi-scalar theory constructed in 91 1 as a contrasting 
ahernative to GR. It depends on two free parameters characterizing how the two scalar fields are coupled 
to matter. Even for strong couplings, the model was fine tuned to pass anyway solar-system tests . Th e 



combination of various binary-pulsar tests however severely constrains both of the parameters 9l|, Il31j . 



44 



hand, like in modified gravity, the matter action will be assumed to depend only on the 
matter fields and their first derivatives, or, in the case of test particles, on their positions 
and velocities only. Moreover, in order to preserve the weak equivalence principle, we will 
assume that all matter fields are coupled to a single second rank symmetric tensor g^j^^. In 
other words, the action of the theory reads 

S = -—^ / V^R* + Sm^tter[lp]g,iu], (5.1) 

iOTTGr J C 

where as before ip denotes globally any matter field. Therefore, we do not seem to have 
changed anything with respect to GR, and this looks like a model in which neither inertia 
nor gravity are modified! However, the crucial difference with GR is that we will assume 
the physical metric Qf^i, to depend not only on the Einstein (spin-2) metric 5^*^, but also 
on the curvature tensors one may construct from it and their covariant derivatives. Several 
definitions of ^^1/ yield the right MOND phenomenology at leading order, and differ only at 
the post-Newtonian level, i.e., for relativistic corrections of order (9(l/c^) with respect to 
the Newtonian and MOND accelerations. Let us choose a particular example: 

~ ^ . , ./^{d,GBfhiX)g;, + 2d,GBd.GB 
^?M^-V+ 4c (0*06/10)7/4 ' 

where GB = R*^up„ — ^R*^^ + Rl is the Gauss-Bonnet topological invariant (which may be 
replaced by -R*^^^ or the square of the Weyl tensor C*^^^. without changing the results at 
first order in a/oq), where □* denotes the d'Alembertian operator (covariant with respect to 
the Einstein metric (?*;,), and where 



The fixed length scale £0 is introduced to make this quantity dimensionless. The numerical 
factors entering Eqs. fl5.2l) - fl5.3p have been chosen to simplify the expression of the function 
h{X), which will be specified below. Note that g^i, is indeed a second-rank symmetric tensor, 
since it is constructed only from tensors, their covariant derivatives and their contractions. 
The matter action S'matteriV'; Qfiu] thereby satisfies exactly Lorentz invariance. Moreover, this 
is still a local action since the physical metric (15. 2p involves a finite number of derivatives 
of It should be noted that Eq. (15.21) does not make sense in fiat spacetime, because it 
would involve ratios of vanishing expressions, but our Universe is never strictly fiat since 
massive bodies do exist. 

Since in Eq. (15. ip . the physical metric Qfj^^ only enters the matter action, the curvature ten- 
sors it involves do not contribute to the dynamics of gravity in vacuum. This implies notably 
that the Birkhoff theorem is satisfied, like in GR. The metric generated by a spherically sym- 
metric distribution of matter is thus of the Schwarzschild type, i.e., reads in Schwarzschild 
coordinates 

glJx^dx" = - (1 - c^dt' + (1 - ' dr^ + (^^' + sin' 6 , (5.4) 

where M is a constant, very close to the one that GR would have defined (but slightly 
modified because of the nonminimal coupling of matter to curvature). It is then straightfor- 
ward to compute the gravitational potential V = — 1(1 + 900)0^ felt by a test particle slowly 
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moving in such a spacetime. The Schwarzschild solution (15 ■4p imphes 



GB = Rl^p„ = C*^p^ = 48 \^-^ j , (5.5) 

□■CB^...o(|^)=(:-H-), 

SO that the time-time component of the physical metric reads 



2GM JGMao , / r \ _ / 1 \ , , 

One thus recovers the MOND potential if one chooses h{X) = InX, and the length scale io 
is then arbitrary since it does not enter any observable prediction. However, we know that 
the MOND potential should not be felt at sma ll dis tances, because it would be inconsistent 



with solar-system tests of Kepler's third law |113l |. One should thus choose a more fine- 
tuned function h{X), such that its derivative gives a 1/r force at large distances r 3> £o but 
a vanishing one when r -C io- One may choose for instance h'{X) = X/{1 + X)^, i.e., after 
integration, 

/i(X) = (l + X)-i+ln(l + X). (5.8) 

Now the constant io does enter the physical predictions, and in order to obtain the correct 
phenomenology, it should be chosen much larger than the solar-system size but small with 
respect to the radius of any galaxy. This underlines that the present model is also quite fine 
tuned, in a similar way as the RAQUAL models discussed in Sees. Ill Fl andl lll Dl above. They 
needed specific aquadratic kinetic terms for the scalar field, and the transition between the 
Newtonian and MOND regimes then occurred at a radius ~ o? \/ GMj a^^ generically much 
too small to be consistent with solar-system tests. Here the transition scale io is fixed by 
hand in Eq. (15.31) . which may be considered as a bonus or a drawback. Note anyway that it 
would be easy to impose io = \/GM/ao, for instance, merely by using again Eqs. (15. 5p and 
(15. 6p to define this quantity and replace it in (15.30 . However, the corresponding definition 
of X ^ (□*GB)5/^v4^/(5aGB)2c would be even more unnatural! 

The light deflection can also be derived straightforwardly from the physical metric (15. 2p 
expanded at first order in Since photons do not feel any global factor of the metric, we 
may actually discard the term proportional to ^/aQg*^J^ in Eq. (15.20 . modulo higher corrections 
O^ao/c^). Let us denote as g^u the resulting metric, which suffices to derive the light 
path. The only remaining modification with respect to general relativity is the final term of 
Eq. (15. 2p . proportional to 9^GB9,^GB. In the case of the Schwarzschild solution (15. 4p . this 
term only contributes to the radial-radial component of the metric 

^ 2GM A^/GM^o ,^(^\ r.nx 

We thus recover the disformal metric constructed in Sec. IIII CI above, and the corresponding 
light defiection, Eq. (13.241) . reproduces exactly (by construction) the result predicted by GR 
in presence of dark matter. Note that the extra (constant) defiection angle li^^GMaoj (? 
is also predicted in the solar system. It would be possible to suppress it by introducing 
another fine-tuned function 6(X) multiplying the c}^GBc}j,GB term in Eq. (15.20 . However, 
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this is not necessary, since this extra deflection angle is ~ 2 x 10~^ arcsec for M = Mq, i.e., 
neghgible with respect to present and foreseen experimental accuracy. 

Therefore, the above model reproduces both the MOND dynamics and the correct light 
deflection in a particularly simple way. It does not introduce any superfluous degree of 
freedom in the vacuum, outside matter sources: It has the same spectrum as GR, namely a 
massless spin-2 fleld carrying positive energy, and does not contain any ghost nor tachyon. 
Moreover, it satisfles general covariance, Lorentz invariance (no preferred frame), locality, 
and the weak equivalence principle. Binary-pulsar tests are also obviously passed, since the 
metric (7*^ generated by a binary system takes exactly the same form as in GR. The only 
difference is that the masses slightly differ from those GR would have predicted for a given 
amount of matter. However, since these masses are not directly measured, but determined 
from the combination of various orbital observables, the analysis of binary-pulsar data is 
actually strictly the same as in GR. 

It should also be underlined that nonminimal couplings to curvature, such as those in- 
troduced in (15. 2p . do appear i n scalar-t ensor theories and in GR itself when describing the 



motion of extended bodies |94l . 11341 . Il35l |. Indeed, when tidal effects can be neglected (i.e., up 
to order included, in GR), compact bodies are characterized only by their mass besides 
their position and velocity, and their action can thus be written Spp = — f mc^ —g*^j^dx^^dx^ . 
As soon as tidal deformations start influencing their motion (from order l/c^° in GR), they 
must be described by more "form factors" characterizing how they are coupled to derivatives 
of the metric. Phenomenological couplings to -Rj^^pg. and higher derivatives are thus expected 
to occur even within GR. Therefore, those entering Eq. (15.21) are rather natural, and deflne 
a consistent dynamics for massive point particles. 

However, the above model also presents several difficulties. First of all, it has been flne- 
tuned to yield the right phenomenology, but other functions of the curvature tensors would 
obviously give a fully different physics. For instance, it suffices to change the numerical 
coefficient of the disformal contribution d^j^GBdyGB in Eq. (15.21) to change the prediction 
for light bending by any amount. Therefore, this model is unpredictive, and just manages to 
fit some experimental data. Actually, this is probably one of its main interests: It illustrates, 
in a very simple framework, all the hypotheses which are also needed in other models of 
the literature. As discussed at the end of Sec. IIIIDI above, the best present MOND-like 
field theory [so], 31, 33| could also predict a fully different light bending than GR if one had 
chosen different functions of the scalar field in Eq. (13.301) . 

A related problem is that there exist many different realizations of the physical metric 
Qfj^u yielding the same phenomenology at lowest order. Indeed, Eq. (15.21) is actually a mere 
rewriting of the needed metric (I5.7p - (I5.9I) in a covariant way. Any other combination of 
the Riemann curvature tensor and its covariant derivatives, reducing to the same metric 
at order 0{l/c^), is thus equally valid. For instance, the crucial factor {2\/GMaQ/ c^)g*^^ 
may equivalently be written in terms of (□*GB)~'^/^/2*p^^i?*^°"^. Any other nonvanishing 
contraction, hke (V*-R*po.^)^ or R^'^^^R'^'^ ^^R'^'^^^, may also be used to build the appropriate 
MOND potential. [On the other hand, neither the Ricci tensor R*^^, the scalar curvature 
R*, nor any of their covariant derivatives, can be used to construct the physical metric 
g^u, since they vanish in vacuum and thereby cannot be felt by a test particle far from the 
gravitational source.] The physical metric (15. 2p that we chose above is convenient to write 
because it involves very few explicit indices, but it has no special property which distinguishes 
it from other possible choices. This freedom actually teaches us another important lesson: 
There is no reason, for a model reproducing the right phenomenology at lowest order, to 
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remain valid at higher orders. In the present case, there is an infinity of ways to complete it 
at higher post-Newtonian orders, by adding arbitrary combinations of the curvature tensor 



and its covariant derivatives. In the case of the TeVeS model [30|, l31|, l33j, there is no reason 
either to trust its higher post-Newtonian predictions. 

Another related problem is that the variable X defined in Eq. fl5.3l) may become imaginary 
in particular backgrounds where 0*06 < 0. One may of course take its absolute value in 
order to always define X consistently, but it would anyway diverge on the hypersurfaces 
where 0*03 = 0. The full metric (15. 2p is actually ill-defined on such hypersurfaces, since 
our "correcting terms" OI^/oq) may blow up. This illustrates once again that the specific 
choice (15. 2p is not justified by any symmetry principle nor any fundamental theory. This 
is just a possible phenomenological realization of the MOND dynamics at lowest order, but 
not a complete theory remaining consistent in all physical situations. It is always possible 
to refine the model to avoid such pathological situations, as we did in Sec. IIIFI to avoid 
the singularity of RAQUAL models on the hypersurfaces where s = g'^'^ d ^ipd ^ip vanishes. 
However, such refinements do not make the models more fundamental, and there is still no 
reason to believe their predictions in situations very different from clustered matter (i.e., 
solar system and galaxies). 

But the most serious difficulty of the nonminimally-coupled model (I5.1l) - (l5.2p is its consis- 
tency within matter. First of all, the physical metric g^^, should always remain hyperbolic, in 
order to define a well-posed Cauchy problem for matter fields. This is clearly the case if the 
correcting terms 0{^/a^) entering (15.21) remain small with respect to the fiat background 
metric 7]^^ = diag(— 1, 1, 1, 1). In more general situations where they become of order 1, 
higher-order corrections are anyway expected to grow too, and the above phenomenological 
model just stops being predictive. Once again, the lowest-order MOND dynamics does not 
suffice to fully specify the physics in all regimes. The deadly problem arises in the field 
equations for gravity within matter. Indeed, since higher derivatives of (7*^, enter the matter 
action, the full Hamiltonian is not bounded by below, because of Ostrogradski's theorem 
discussed in Sec. Ill CI above. Therefore this model is unstable, in the same way as higher- 
order gravity. Although one cannot identify any ghost degree of freedom around a fiat empty 
background, a Lagrangian of the form R + 7GB + a{Rf^^p^)'^ {n > 2, and probably even 
i?-t-/(GB) as in Refs. 73 , 74 , 7^]) does contain negative-energy modes around curved back- 
grounds, and the present model (I5.1l) - (l5.2p involves negative-energy modes within matter. 
This is illustrated in terms of Feynman diagrams in Fig. O 

As a simple toy model of such unstable theories, one may consider for instance the 
Lagrangian C = —{d^ip)"^ [1 + '^{di,ipY] in fiat spacetime, where A is a nonzero constant. 
The scalar ip is an analogue of the matter fields in our present model (I5.ip - (l5.2p . and of the 
usual massless graviton in Lagrangians of the form i?+7GB+a;(i?^j,^^)". One cannot identify 
any kinetic term for the second scalar (f around a background ip = const, but it is obvious 
that it defines a ghost degree of freedom in backgrounds where Xld^ipy < 0. More generally, 
one can easily compute the Hamiltonian n= [V^^ + [di-^f] [l + X{dj(py] + A^^ [{diipf - 3tp^] , 
and prove that there always exist initial conditions such that Ti. —00. If A > 0, it suffices 
for instance to choose diip = di^p = 0, 0^ = 1/A, ^/i)^ — > 00, and if A < 0, to choose dip = 0, 
= 0, {diipY'p'^ — *• 00. Both cases can be expressed the same way in terms of the conjugate 
momenta vr^ = dC/dip and vr^ = dC/dip: Whatever the sign of A 7^ 0, one gets H, — >• —00 
for = 0, vr^ = 0, diip and vr^ — > 00. 

In the present model (I5.ip - (l5.2p . the simplest way to understand the existence of negative- 
energy modes, confined within matter, is to consider a homogeneous density (either fill- 
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(a) -r^Kmr; = but ^nrfC... ^ 



(b) -KKKKK = but _ 

FIG. 6: Diagrammatic illustration of the negative-energy modes spoiling the stability of (a) La- 
grangians of the form R+jGJi + a{R'j^j^p^)^ , n > 2, and (b) the present nonminimal metric coupling 
(j5.ip - (|5.2p . As before, curly lines represent spin-2 degrees of freedom, namely the usual massless 
and positive-energy graviton in plain line, and the massive ghost in dotted lines. Contrary to FigjSl 
straight lines represent here matter fields V instead of the scalar field ip. In both cases (a) and 
(b), no ghost degree of freedom can be identified around a flat or empty background. However, 
a ghost mode does exist in curved backgrounds in case (a), and within matter in case (b). They 
use the background gravitons or matter fields as "catalyzers" in order to propagate and develop 
instabilities. 



ing the Universe, as in standard cosmological models, or in the interior of a massive 
body). Then the matter action Sraa.tterl^'', 9 fj.u] actually defines a gravity theory of the form 
f(R*„,^„, R*,, n ^X R* ,,,f.^^ ■ ■ which is known to contain a massive spin-2 ghost degree of 
freedom 69|, [TIJ, [72], as recalled in Sec. Ill C[ Of course, one may invoke again higher-order 



corrections to the phenomenological model (15. 2p . which may be able to stabilize it. However, 
it remains that it is unstable without considering such higher-order corrections, and there- 
fore that its lowest-order truncation (15.20 cannot describe consistently massive bodies. At 
the classical level, one should analyze carefully this instability to determine its characteristic 
timescale. Naively, one could expect it to be related to the MOND acceleration scale oq, 
i.e., of order c/ao ~ Q/Hq 8 x 10^" yr. Since this is several times larger than the age 
of the Universe, this would obviously not be a problem for massive stars. However, other 
dimensionful constants do exist in the model. The length scale C.q entering (15. 3p does not 
change any physical prediction at small enough distances, but within matter, all masses 
or coupling constants entering the matter action S'^atter may contribute to the instability 
timescale. Therefore, the analysis of this instability is highly nontrivial. It may actually 
happen that it does not develop at the classical level. Indeed, it is quite difficult to exhibit 
toy models in which the presence of ghost degrees of freedom (negative kinetic energy) im- 
ply classical instabilities, whereas they are obvious in presence of tachyons (maximum of 
a potential). On the other hand, ghost degrees of freedom are known to be deadly at the 
quantum level, because the vacuum then disintegrates instantaneously (70| . 

In conclusion, the nonminimally-coupled model (I5.1l) - (l5.2p has many interesting features, 
notably the fact that it is rather simple and that its predictions can be easily computed 
from the general relativistic metric. But it is unstable and therefore one cannot trust its 
predictions. It may still be useful as a phenomenological framework at lowest order in 
powers of ^/oq, since it goes beyond the original MOND formulation [13]. For instance, 
when a test particle is located between two massive bodies, it may happen that its total 
acceleration is smaller than oq although each individual force happens to be in the standard 
Newtonian regime (/ > mao). In such a case, one may wonder what MOND should predict, 
whereas the present nonminimally-coupled model unambiguously defines the dynamics from 
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the local curvature tensor and its covariant derivatives. This model also allows us to compute 
consistent MOND effects for non-spherically symmetric sources. But its main interest is 
to show that the MOND dynamics is consistent with most of GR's symmetries. It also 
illustrates in a clear way many generic problems of other MOND-like field theories of the 
literature, notably their unpredictiveness (because they are just tuned to reproduce the 
needed physics) and their almost systematic instability. 



B. A nonminimal scalar-tensor model 

The most serious problem of the above nonminimally-coupled model is its instability be- 
cause of the presence of higher derivatives of the metric in the matter action. By introducing 
back a single scalar field, we will now show that this instability may be avoided (but that 
a more subtle consistency problem anyway remains, as discussed at the end of the present 
section). It will not even be necessary to give it an aquadratic kinetic term as in Sec. IIIEI 
A standard scalar-tensor theory of the form (12.121) will suffice for our purpose, but with a 
more complex physical metric g^j^iy. 

The idea is very similar to the previous section IV^ We showed in Eqs. (I5.5p - (I5.6I) that a 
combination of curvature tensors and their covariant derivatives gave us a local access to the 
baryonic mass M and to the distance r, independently from each other. It was thus possible 
to construct from them the needed MOND potential y/GMao In r. In standard scalar-tensor 
theories, with an exponential matter-scalar coupling function A{(p) = e°"^, we know that the 
scalar field takes the form ip = —aGM/rc^ + 0{l/c^) in the vicinity of a spherical body of 
mass M, if the background value (fo vanishes. Then it is obvious that combinations of ip 
and its derivative dr^p = aGM/r'^c^ + 0{l/c^) also give us access to M and r independently. 
One may thus imitate the previous model (I5.1l) -( l5l2i) by defining 

f d'^x . \ R* s 1 /mc\2 A r, ~ n .^n 

S = ^j — \T - 2 " 2 ^ / + ' ^^-^^^ 

where s = g^^d^ipd^ip denotes as before the standard kinetic term of the scalar field, and 
where g^y = A'^g^^ + Bd^ipd^ip takes the disformal form (13.211) with 

A((^,s) = e°'^ - ^/i(X), (5.11a) 
a 

I . , ^ 

Bi^.s) = -A- , (5.11b) 

a s 

X ^ v^,-V4. (5.11c) 

c 

The function h{X) = (1 -|- X)~^ + ln(l -|- X) was chosen in Eq. (15. 8p to reproduce the 
MOND 1/r force at large distances but a negligible one at small distances. The potential 
V{(p) = ^{mc/h)^(p'^ introduced in Eq. (I5.10p has a negligible influence on local physics 
(solar system, binary pulsars, galaxies and clusters) if the mass it involves is small enough, 
say mc^/h ^ H where H denotes the Hubble constant. However, it plays a role in the 
cosmological evolution of the scalar field, driving it to the minimum ipo = 0. This term 
ensures thus that the background value of the scalar field vanishes. Note however that the 
above model must be refined to remain consistent in the cosmological regime (s < 0), on the 
hypersurfaces s = surrounding any galaxy or cluster (see Sec. Ill Fl above), and near the 
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center of any massive body where s — 0. It suffices, for instance, to replace s by 
in Eqs. flS.lip . where e ~ {mc/flY ^ H^/c^ is small enough to have negligible influence 
both on Newtonian and MONDian physics. As already underlined below Eq. (12.181) and in 
Sec. IV At such a refinement just illustrates how to define a consistent field theory, but there 
is no reason to trust its predictions in the cosmological regime s < 0. Indeed, many other 
functions of s could connect smoothly to those entering (15.111) for small values of s, and they 
are not imposed by the MOND phenomenology. 

Since the above model defines a standard scalar-tensor theory outside matter (in vacuum), 
we know that the scalar field takes the form Lp oc —GM/rc^ + 0{l/c^) in the vicinity of 
a spherical body of mass M. The multiplicative factor is just a in scalar-tensor theories, 
but it may now be modified by the MOND-like corrections we added in the physical metric 
cjfj^iy. One expects it to remain close to a if the matter sources are compact enough to be 
in the Newtonian regime of the model. This is not only the case for solar-system bodies 
and binary pulsars, but also for almost all the baryonic mass of galaxies and clusters. We 
will thus merely assume that ^ —aGM/rc^ in the present discussion. We will come back 
below to the (crucial) problem of the scalar field equation within matter, and thereby to the 
actual coupling constant a entering this expression. The value cp ~ —aGM/ rc^ implies that 
X = r/ ^/GMJoq, which not only gives us a direct access to the radial distance r, but also 
involves the MOND characteristic radius £o = ^/GMJoq. This is thus already a bonus with 
respect to the purely metric model (I5.1l) - (l5.2p . since one needed to impose by hand that the 
arbitrary scale io was much larger than the solar-system size but smaller than the radius 
of any galaxy. Here, this is obtained automatically from a rather simple definition of X, 
Eq. fISTTcD . One also finds that ipX/a = ~y/GMa^/c^, i.e., precisely (by construction) the 
needed factor of the MOND logarithmic potential. Therefore, the physical metric admits 
the following post-Newtonian expansion: 



^ 2a^GM 2VGMao , 
1 — + 2 ^ 



A^GM^o 



6;6: + o 



1 



- . (5.12) 



For distances r ^ a/ GM/oq, the function h{X) has been chosen to reproduce the MOND 
logarithmic potential, and the disformal contribution oc 5^5^ in (I5.12p is precisely the one 
needed to reproduce the same light deflection as GR in presence of dark matter, Eq. (13.241) . 
as already noticed below Eq. (15.90 for the previous nonminimal metric model. On the 
other hand, for r <^ a/ GM/qq, the function h{X) tends to 1 and does not contribute to 
any physical observable, whereas the disformal contribution oc 5^5^ still gives the constant 
anomalous light deflection (I3.24p . but too small to be of observational signiflcance in the 
solar system. Therefore, this model does reproduce the right phenomenology. Solar-system 
tests just impose the tight bound < 10~^, as in standard scalar-tensor theories, because 
the —2Q^GM/rc^ term contributes to our interpretation of light deflection, as explained in 
Eq. ([32D above. 

Since this model has been devised to reproduce the standard Brans-Dicke predictions 
for distances r <^ GM/aQ, we can al so conclu c le th at it passes binary-pulsar tests, which 



impose the constraint < 4 x 10 ^ 129 . ISOl . 131 |. weaker than the above solar-system 



This negative value ensures that the matter field equations are well posed, since inequalities p.22p are 
satisfied: + sB = + Ay/GM^/c^ > > 0. 
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one. Contrary to our discussion of the TeVeS model in Sec. IIVDI the (spin-2 and spin- 
0) gravitational waves are now strictly the same as in Brans-Dicke theory, even at large 
distances where MOND effects can be observed on material test bodies. Indeed, this model 
is a pure scalar-tensor theory in vacuum, so that gravitational waves themselves do not 
feel any MONDian effect. This is an interesting feature of this model, and also its main 
difference with the previous ones proposed in the literature. In previous models, the helicity- 
(scalar) waves undergo MONDian effects at large distances from massive sources, whereas 
they follow the same geodesies (of the Einstein metric g*^^) as the helicity-2 waves in the 
present model. Unfortunately, this different prediction cannot be used to experimentally 
discriminate between these models, because the solar-system constraint < 10~^ on the 
matter-scalar coupling constant implies that gravitational-wave detectors (made of matter) 
are almost insensitive to scalar waves. On the other hand, let us underline that they could 
in principle discriminate between such MOND-like field theories and the standard dark 
matter hypothesis. Indeed, all RAQUAL and disformal models, including TeVeS and the 
present one, predict that the helicity-2 waves are not directly coupled to the scalar field, and 
thereby are not deflected by what we usually interpret as dark matter haloes, contrary to 
light (coupled to the physical metric c/ni, describing MONDian effects). Therefore, when the 
LISA interferometer will detect gravitational waves emitted from optically known sources, 
their direction should be slightly different from the optical direction if there exists at least 
one intermediate gravitational deflector. Unfortunately again, the LISA angular resolution 
will be about the angular size of the Moon (30 arcmin), much too large to measure such 
a deflection angle between electromagnetic and gravitational waves. But the time- delay 
between gravitational waves and optical or neutrino pulses from supernovae would all ow u s 



to discriminate between MOND and the dark matter hypothesis, as recently shown in [136 



The deadly problem of the previous model fl5.1l) - fl5.2p was its Ostrogradskian instabil- 
ity. One may thus wonder if the present scalar-tensor variant avoids it. Indeed, since the 
matter action a priori involves covariant derivatives of the matter fields, it depends on the 
first derivatives of the physical metric c/f^i,, and thereby of the second derivatives of (f, see 
Eqs. (15. lip . The same instability as before, caused by second and higher derivatives, seems 
thus to spoil the present model too. However, let us recall that the standard model of 
particle physics only involves gauge bosons described by 1-forms, and fermions described by 
Dirac spinors. If is well known that all covariant derivatives may be replaced by ordinary 
ones in forms, because their antisymmetry cancels all the Christoffel symbols T^^. There- 
fore, there is actually no derivative of Qf^u involved in the action describing the dynamics 
of gauge bosons, and only ip and its first gradient dfj,ip enter it. On the other hand, the 
action of a spinor in curved spacetime does dep end of the first derivative of the metric (via 
the derivative of a tetrad, cf. Sec. 12.5 of |l22 |). but only linearly. Schematically, it reads 
■il)g'^(^ip + ipg^{(^g)%lj. Therefore, V*^dj^(p enters this action, but if one tries to define a conju- 
gate momentum p2 = dC/dip, then its expression does not depend on any second derivative 
of if. Therefore, one cannot express </3 as a function of (f, (p and p2 as in Eq. (12.41) . and we are 
thus precisely in the case of a degenerate Lagrangian for which the Ostrogradskian instabil- 
ity does not occur. More explicitly, all second time derivatives of the scalar field, (p, may be 
eliminated by partial integration. In conclusion, the action of the standard model of particle 
physics actually does not depend on (p, and the above model is thus free of Ostrogradskian 
instability. 

However, this absence of generic instability does not suffice to prove that the model is 
indeed stable. One must also check that the scalar field equation remains consistent within 
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matter. In the most general case of arbitrary matter fields and arbitrary initial conditions, 
this problem is quite involved because the matter action depends on the derivative of the 
scalar field. As discussed in Sec. IIII Ct one should write all the field equations, including 
for matter, and prove that the second derivatives of the scalar field define a hyperbolic 
operator. The simpler case of a pressureless perfect fluid will however suffice to exhibit the 
subtle difficulties appearing in this class of models. The generalized RAQUAL kinetic term 



written here in the rest frame of the perfect fluid and in coordinates locally diagonalizing the 
Einstein metric (7*^ = diag(— 1, 1, 1, 1), should satisfy inequalities (al) and (bl) quoted below 
Eq. (13.271) . The contribution of the scalar field to the Hamiltonian should also be bounded 
by below. Since the model has been constructed to reproduce the MOND dynamics at the 
lowest (Newtonian) order, without imposing its behavior at higher post-Newtonian orders 
nor in the cosmological regime, let us focus on the conditions (a2), (b2) and (c2) that / 
should at least satisfy at the Newtonian level. As mentioned in Sec. IIII C\ condition (c2) is 
actually implied by (a2), because we imposed a positive value of the coupling function B in 
order to increase light deflection in the MOND regime. Moreover, one can check easily that 
the terms proportional to (c^qV?)^ may be neglected in (a2) and (b2). On the other hand, 
note that A'/c^ and sA" /c^ are actually of Newtonian order (9(l/c°), in spite of the explicit 
factor At this Newtonian order, it suffices thus to check conditions 



where / may be approximated as / ~ s + {87rGp/c^)A, as if we were considering a RAQUAL 
model in vacuum; see conditions (a) and (b) of Sec. Ill El In such a case, these hyperbolicity 
conditions suffice to ensure that the scalar-field contribution to the Hamiltonian is bounded 
by below. [The coupling function B plays an even more negligible role in the Hamiltonian 
than in the hyperbolicity conditions, namely of order 0{l/c^) smaller than the lowest term 
involving the matter density p.] 

In the Newtonian regime, i.e., when the MONDian correction involving ipX/a in 
Eqs. (15. lip are negligible, it is obvious that the hyperbolicity conditions (a2) and (b2) 
are satisfied. The model reduces indeed to the standard Brans-Dicke theory, whose consis- 
tency is well known. If all matter in the Universe was clustered as compact enough bodies, 
like stars and planets, then this Newtonian regime of the model would always be reached 
when considering the scalar field equation within matter, and its Cauchy problem would 
always be well posed. As mentioned below Eqs. (15. lip , there would still remain a difficulty 
near the center of any body, where the spatial gradient of the scalar field tends to 0, and 
therefore where MOND corrections may become large (before they are saturated by our ad 
hoc replacement of s by a/s^ + e"^) . However, it is not even necessary to study such centers 
of compact bodies to point out an extremely serious difficulty of this model. Indeed, there 
does exist dilute gas in the outer regions of galaxies where the MOND phenomenology ap- 
plies. One should therefore check whether conditions (a2) and (b2) hold when the (pX/a 
corrections start dominating even over the Newtonian force. Actually, for the coupling func- 



tions (15. lip , one finds that condition (b2) is easily satisfied, but it happens that (a2) is 
not. Indeed, /' ^ 1 + (STrGp/c^)^', where A = {h + X dh/dX) (fX/Aas. The coefficient 




(5.13) 



(a2) f > 0, 



(b2) 2sf" + /' > 0, 
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{h + X dh/dX) varies from 1 to 3.3 in the range of distances relevant to the MOND regime, 
namely X = r/^GM/ao < 10. On the other hand, ^X/as = -{^/GMao/c^) {aGM/r^'^)-^ 
is a huge negative number, because the kinetic term s of the scalar field appears in its de- 
nominator. Even multiplied by the small gas density Gp/c^ existing in the outer regions of 
galaxies, this negative contribution to /' dominates by several orders of magnitude over the 
1 coming from the standard kinetic term in vacuum. In other words, /' < within such a 
gas, and the scalar perturbations are thus ill defined. In spite of its nice features and its 
relative simplicity, the model (15.111) is thus inconsistent, because the scalar field equation is 
not always hyperbolic within matter. Such a difficulty was not discussed in previous works 
on MOND-like field theories, although it seems generic for disformally-coupled scalar-tensor 
models (fMH). 

Of course, one may try to tune the coupling function A, Eq. (15.1 lap , to cure this incon- 
sistency. For instance, since (pX/a takes the constant value —y/GMao/c^ outside matter, 
one may add to A any function of ipX/a without changing its spatial gradient, i.e., the force 
felt by a test mass. However, it is straightforward to show that there is not much freedom 
for such a function. In order to ensure condition (a2) without spoiling (b2), it must be close 
to k{ipX / a)""^ , where /c is a constant. This precise expression has the peculiarity not to con- 
tribute to condition (b2), whereas it adds a positive contribution to /'. The crucial difficulty 
is that the proportionality constant k must take a fixed value, imposed in the action of the 
model, and should not be chosen independently for each galaxy or cluster (cf. our discussion 
in Sec. Ill Bl above). To guarantee condition (a2) in any situation, k should thus be chosen 
large enough, and one finds that k > 2(G'Mciusterao/c^)^'^^ would suffice, where Mduster is the 
largest cluster mass existing in the Universe. Obviously, this appearance of another dimen- 
sionful constant in the model underlines that it would be even more fine tuned than before. 
[One might be tempted to simply set = 1 to make the model more natural, but the coupling 
constant A would then take large values instead of admitting a post-Newtonian expansion 
A = 1 + (9(l/c^), and the predicted MONDian effects would then be reduced by the inverse 
of this large value.] But it happens that such a tuning of the coupling function A actually 
does not work. Indeed, although a term klipX/a)'"^ would not contribute to the force felt 
by a test mass, it does contribute to the scalar field equation 2V*(/'V^V9) = df /dip. In 
the MOND regime, it even dominates over the other contributions in both sides of this field 
equation, in spite of its coefficient (SnGp/c^) involving the small gas density p. In a region 
where p is almost constant, one may thus divide each side by it, and one gets a scalar field 
equation fully different from the one satisfied in vacuum,^^ although it does not depend any 
longer on the precise value of p. Its exact solution is difficult to find, but the signs involved 
anyway show that (p would quickly tend to within matter (in the MOND regime), instead 
of keeping the form ip ^ —aGM/ rc^ that we assumed from the beginning of this section to 
construct the present model. Even without invoking the signs on the different terms in this 
field equation, it is anyway obvious to check that cp ~ —aGM/rc^ cannot be valid within 



This problem actually already arises with the initial coupling function A, Eq. (|5.11ap . whose MONDian 
terms dominate within the gas in outer regions of galaxies. However, one can check that the scalar field and 
its derivative keep the same orders of magnitude as the vacuum (Brans-Dickc) behavior ip « ~aGM/rc^ 
we assumed, so that the phenomcnological predictions remain correct up to factors of order unity. On 
the other hand, the scalar field behavior is fully changed within this outer gas when the correcting term 
k{(pX / a)~'^ is included. 
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the gas in outer regions of a galaxy. In other words, the extra force predicted by such a 
model, if any, would not take the MOND form that we wished. In conclusion, a tuning 
of the coupling function A to ensure the hyperbolicity of the scalar field equation within 
matter suffices to ruin the predictions of the model. 

The above hyperbolicity problem underlines again how difficult it is to construct a 
MOND-like field theory, even fine tuned. In spite of the presence of several free functions 
in the model, which could a priori fit any phenomenology one wishes, the consistency of 
the field equations within matter is a strong enough condition to rule it out. One should 
however underline that this conclusion cannot be considered as a no-go theorem. The same 
framework as Eqs. fl5.10l) - fl5.11l) might be able to reproduce consistently the MOND dynam- 
ics with different coupling functions A and B. If the scalar field does not take the form 
—aGMjr(? in the MOND regime, then specific combinations of (/J and its kinetic term s 
might be able to give us access independently to the crucial factor \J GMao / (? and to the 
radius r / a/ GM/aQ. The MOND metric (15.121) might thus be obtained again, and our dis- 
cussion above shows that the Ostrogradskian instability would still be avoided. We expect 
the hyperbolicity conditions (a2) and (b2) to still cause serious problems within matter 
in the MOND regime, but we do not have any proof. Our investigations indicate anyway 
that the coupling functions A and B would very probably involve unnaturally complicated 
expressions (like large powers of s or transcendental functions). 



C. The Pioneer anomaly 

As mentioned in Sec. HI the two Pioneer spacecrafts exhibited an anoma lous extra accel- 
eration 5a ~ 8.5 X 10~^°m.s~^ towards the Sun, between 30 and 70 AU [l^, 137 . 138 |. From 



a theoretical point of view, this anomaly is a priori much easier to explain than dark matter 
by a fine-tuned model modifying Newtonian gravity. Indeed, it has been observed in our 
single solar system, mainly characterized by the solar mass M©, whereas galaxy and cluster 
rotation curves have confirmed the relation Mdark oc Mbaryon for many different values 
of the masses. A possible explanation (or, rather, description) of the Pioneer acceleration 
would be to couple matter to a scalar field, whose potential is merely fitted to reproduce 
observed data. However, in order not to spoil the precision tests of Kepler's third law in the 
solar system, such a potential should not manifest at distances smaller than Saturn's orbit 
{k, 10 AU), and cause an almost constant extra force beyond 30 AU. Besides the fact that 
this would be quite unnatural, one may also wonder if this would correspond to a stable 
theory. Indeed, the simplest stable potential for a scalar field is a mass term cx <y9^, and 
it causes a force which decreases (exponentially) as the distance grows. Even within the 
a priori simpler case of the Pioneer anomaly, the construction of a stable model does not 
seem obvious, thus. 

A first, important, step has been understood in the series of papers [l39l . [l40l . [mi] 



by Jaekel and Reynaud. From a phenomenological viewpoint, it is much more natural to 
describe the Pioneer anomaly by a modification of the spatial metric Qij rather than the 
Newtonian potential involved in qqq. Indeed, if the metric is written in the Schwarzschild 
form (13.121) . then the geodesies equation for a test mass mo reads 

S^/c^ - {mlc^ + J\^) B = J^ABu'\ (5.14) 

where £ is the particle conserved energy (including rest mass moc^), J its angular mo- 
mentum, u = 1/r, and a prime denotes here derivation with respect to the angle in polar 
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coordinates. For a circular orbit, m' = so that the spatial metric component grr = A 
does not enter this equation. More generally, its right-hand side is proportional to the 
square of the orbit's eccentricity, e^. In conclusion, planets on nearly circular orbits are 
almost insensitive to a small modification of A, whereas they would feel any extra potential 
entering the time component B. On the other hand, hyperbolic trajectories,^^ like those of 
light rays and the Pioneer spacecrafts, are directly sensitive to a modification of A. This is 
actually the crucial difference between galaxy rotation curves (which exhibit a discrepancy 
with Newton's law for the circular orbits of outer stars) with the Pioneer anomaly (occur- 
ring for hyperbolic trajectories). While it was necessary to modify the time-time metric (700 
to reproduce the MOND phenomenology (together with the spatial one gij to predict the 
correct light deflection), it is possible to account for the Pioneer anomaly by a modiflcation 
of Qrr alone. 

However, post-Newtonian tests in the solar system, notably those of light deflection and 
Shapiro time delay, severely constrain the magnitude of the deviations from general relativity 
in the spatial component A. In particular, it is impossib le to account for the Pioneer a nomaly 
within the parametrized post-Newtonian formalism [l|, 142, 143 . 144, 145 . 146l . 147l |. which 
describes all metrically-coupled theories of gravity at order 0{l/(?) with respect to the 
Newtonian force, under the hypothesis that the gravitational fleld does not involve any 
characteristic length scale. This hypothesis implies that all post-Newtonian parameters, 
including the famous Eddington parameters [3 and 7, must be independent of the radial 
distance r. Within this wide class of theories, the Pioneer anomaly is merely inconsistent 
with other precision t ests i n the sol ar system. 

The lesson of Refs. jl39l . 140 . 141 1 is that a spatial dependence of the Eddington parameter 
7, going thus beyond the PPN formalism, suffices to account for the Pioneer anomaly without 
spoiling the classic tests. This parameter should of course remain very close to the general 
relativistic value 7 = 1 at small radii, so that light-defiection and time-delay predictions are 
consistent with observatio n, but it may change at distances ^ 30 AU to fit Pioneer data. In 
their first analysis, Refs. 139 . 14Cll | showed that an expression 7 = 1 + /tr^ would a priori 
suffice, for an appropriate constant k tuned to reproduce the Pioneer extra acceleration. In 
other words, the radial component of the metric should take the form g^r = l+2'jGMQ/rc^ + 
(9(l/c^) = 1 + (2GM0/c^)(l/r + K,r) + (9(l/c^), in which a small term linear in r corrects 
the general relativistic 1/r potential. The refined discussion of Ref. [l4:l] actually leads 
to a more complex radial dependence of the parameter 7, cubic in r, but the qualitative 
conclusion remains the same. 



Let us just summarize here the simpler result of [139|, Il40l |. The extra acceleration 6a. 



Beware that in the Schwarzschild coordinates we are using, the post-Newtonian term in the time compo- 
nent goo of the metric is proportional to the combination /3 — 7 of the Eddington parameters, whereas it 
is proportional to 7 in the radial component. It would thus be misleading to claim that circular orbits are 
merely independent of 7. The correct statement should underline that this is only true for a fixed value 
of /3-7. 

Recall that the virial theorem implies w^/c^ = 0{GM /rc^) for bound orbits, whereas hyperbolic ones 
reach a nonzero asymptotic velocity jc? 3> GMjrc? for large distances r. As compared to the time- 
time contribution gQQC^dt^ to the line element ds^, the spatial one gijdx''dx^ = gijV^v^dt^ is thus of 
post-Newtonian order for bound orbits, but much larger for hyperbolic ones (and actually of the same 
order as the time-time contribution for hght rays, yielding the famous (I + 7) factor of the hght deflection 
formula) . 
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caused by a term fc„r" added to grr = ^ is a straightforward consequence of the geodesies 
equation (15.141) . One finds 

6a.= --knur"-'^ v"^ n, (5.15) 

where v denotes the particle's velocity and n = x/r the radial unit vector (pointing away 
from the Sun). Since the velocities of the Pioneer spacecrafts were almost constant on the 
range of distances where their anomalous acceleration was observed, one thus needs n = 1 
and a positive value of ki to obtain a constant 6sl directed towards the Sun. Actually, one 
does not directly observe the acceleration of the spacecrafts, but deduce it from Doppler 
tracking. One should thus also take into account the fact that electro magnetic w aves travel 



in the same metric, and feel the anomalous potential too. References |139l . |140| show that 
this leads to a factor 2 difference in the prediction of the recorded anomalous acceleration, 
i.e., Sa = kiv"^. Numerically, for the Pioneer spacecrafts, v ~ 1.2 x lO^m.s"^ and 5a ~ 
8.5 X 10~^°m.s~^, therefore 

k, ^ ^ (lO^AU)-^ (5.16) 

If interpreted as a varying Eddington paramet er 7 = 1 + /tr^, with IkGMqI (? = ki, this 



gives K !^ 45 AU ^. The refined analysis of Ref. |141[ | yields a quadratic extra potential k2r 



,2 



with ^2 ~ 8 X 10~^ AU~^, i.e., a varying Eddington parameter 7 1 + 4(r/AU)'^. It further 
tunes this extra potential in the form k2{r — r^)^ + k2{r — r^), where is the radius of the 
Earth's orbit and fcg another constant characterizing the anomalous acceleration at r = r®. 
We shall not take into account such refinements in the following, but just mention that they 
would not change our conclusion. 

The last difficulty is to exhibit a consistent field theory reproducing such a radial depen- 
dence of 7, both stable and admitting a well-posed Cauchy problem. Because of the generic 
difficulties with modified gravity models recalled in Sec. [Tll and notably the stability issues 
of higher-order g ravity dis cussed in Sec. Ill CI it is a priori not obvious that such a theory 
exists, and Refs. 139l . 140 1 actually write nonlocal field equations without discussing these 



crucial mathematical requirements. But it happens that the scalar-tensor framework of 
Sec. IV Bl above suffices to reproduce consistently a varying Eddington parameter 7. Indeed, 
the disformal contribution Bd^ipd^ip to the physical metric (13.211) is precisely what is needed 
to change only the radial component grr, i.e., what is denoted as A in the geodesies equation 
(I5.14p . As in Sec. IV Bt one may thus consider a scalar field weakly coupled to matter via a 
Brans-Dicke conformal factor A'^{ip) = exp{2aip), and build the appropriate correction^^ to 
the radial metric by using ip ^ —aGMQ/rc? and ^/s = drip ~ aGMQ/r'^c^ . Since we do not 
have any experimental evidence that the Pioneer anomalous acceleration is related to the 
solar mass (nor its square root nor any function of it), several expressions can actually be 
used to build an extra gravitational potential linear in r, for instance 1/cp or s~^/^, or any 
function like ip~^ f {ip"^ / s) . The simplest choice would be to impose B = —X^/ip^, where A5 
is a positive constant, i.e., a physical metric of the form 

~g,u = e''^^g;,-^d,ipd,ip, (5.17) 



A correction to the time-tinie metric ^00 would also be easy to define by using ip and dr^. This might be 
necessary to avoid too large modifications of the planets' perihelion shifts, which depend on both of the 
PPN parameters /3 and 7. 
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so that grr = e^^'^f?*^ — {\5/(p){s/ip^) involves a pos itive correction proportional to r. [The 



more complex, quadratic, potential derived in Ref. |14l[ | can be obtained in the same way. 



by choosing for instance B = fi^/ip^ + fi^/ip^ + fi4/(p^, where /i4,5,6 are three constants such 
that 5 > 0, in order to predict an extra acceleration directed towards the Sun.] Of course, 
expression fl5.17l) is valid only in the vicinity of the Sun, where s > and (f < 0. It should 
be refined to remain valid in other regimes, but the Pioneer data do not tell us how. Since 
the condition + sB > must be satisfied to ensure the hyperbolicity of the matter field 
equations, see Eq. fl3.22p . it might be better to consider coupling functions B = Xgs/\ip\^ 
or B = Xi/{\ip\s). However, it would also suffice to multiply the above B = —X^/ip^ by 
a function h{s), such that h{s) = for s < and h{s) = 1 for s > {aG Mq / r'^^^c^)'^ , 
where rmax ~ 70 AU is the largest distance for recorded Pioneer data. The above met- 
ric (15.171) should also be refined to avoid divergences when ip passes through zero with a 
nonvanishing derivative. A possible solution would be to replace —(p by where 
e <^ aGMQ/rj:aaxc'^ is a tiny dimensionless number, negligible for the Pioneer phenomenol- 
ogy but eliminating our anomalous potential at larger distances. Numerically, one would 
need e < IQ-^^a < 3 x 10"^^ if one takes into account the solar-system constraints on the 
matter-scalar coupling constant a < 3 x 10"'^. On the other hand, to reproduce the observed 
Pioneer anomaly, the constant A5 entering Eq. (15.1 7p should take the numerical value 



A5 ^ a^(lO-^m)^ < (2 X lO^^m)^ (5.18) 

the upper limit corresponding to the largest allowed value for a. [The other possible ex- 
pressions for B quoted above would need factors /ig ~ a^(4 x 10~^m)^ < (4 x 10~^'*m)^, 
Ag ^ a5(0.4m)^ < (3 x lO'^m)^ and Ai ^ lO'^^a < 3 x 10"^^ respectively.] 

Such numerical constants underline that we are merely fine-tuning a model to account 
for experimental data. However, the crucial point that we wish to stress is that the model 
(I5.10p with a physical metric (15.170 is stable and admits a well-posed Cauchy problem. The 
hyperbolicity conditions (al) and (bl) of Sec. IIII CI are indeed trivially satisfied, because 
both coupling functions A = exp{a(p) and B = —X^/ip^ are independent of the kinetic term 
s, and because B > in the regime where this model is defined. Actually, for such a simple 
case, one may even check the hyperbolicity of the fully general scalar field equation (I3.25P 
within any kind of matter. The fact that the Hamiltonian is bounded by below is also an 
obvious consequence of the positivity of B. [Note that all these consistency conditions are 
also obvious for the more complex expression B = fie/^fil + f^s/v^ + fii/'p'^ mentioned above 
to reproduce the quadratic potential derived in Ref. jl4l| . On the other hand, they are 



less trivial for the other possible expressions for B quoted above, such as B = A9s/|<y9|^ or 
B = Ai/(|v9|s). However, they may still be easily checked at the lowest post-Newtonian level, 
where the hyperbolicity conditions reduce to (a2), (b2) and (c2) of Sec. IIII Ct and where all 
terms proportional to u'^d^ip may be neglected.] 

The final consistency check is to show that the scalar field takes indeed the form cp « 
—aGMQ/rc^ in the solar system, since we assumed this expression to build our anomalous 
potential. One just writes the scalar field equation within matter, 2V*(/'V(^v') = df /dip, 

where / is defined in Eq. (I3.27p . If the Sun were an isolated body, the scalar field would 
be generated by it only, and one would have u'^d^ip = strictly (up to a negligible time 
variation of the scalar field possibly imposed by the cosmological evolution). In such a case, 
the model reduces to Brans-Dicke theory, and one gets thus obviously (p ~ —aGMQ/rc? . 
When considering other bodies in the solar system, the largest effects due to the disformal 
coupling function B occur at large distances from the Sun. One finds that the coupling 



58 



constant a is then renormalized into a + 0(tv^ j GMq^kxT j where k\ ~ (lO^AU)"^ is 
the coefficient fl5.16p of the hnear anomalous potential we added to the radial metric. The 
coefficient 0{rv'^ / GMq) is of order unity for planets because of the virial theorem, and solar- 
system tests impose < 10~^. For the largest allowed value of a, we can thus conclude 
that our anomalous potential does not contribute significantly to the scalar field equation 
up to r ~ 10 AU (i.e., Saturn's orbit). For outer planets, a may be increased by a factor 
< 4, but since their masses are negligible with respect to that of the Sun, the approximate 
value of the scalar field near the Pioneer spacecrafts remains anyway ~ —aGMQ/rc? . The 
largest modification of the matter-scalar coupling constant occurs for the Pioneer spacecrafts 
themselves (or any object on a hyperbolic orbit at similar distances from the Sun), since 
rv'^/GM.Q 10 for r = 70 AU and the spacecrafts' velocity ^ 1.2 X 10^ m.s-^ This 
means that this coupling constant is increased by a factor ~ 80 at such distances, explaining 
incidentally why these spacecrafts are sensitive to the scalar field. But the local perturbation 
of ip they induce does not act on their own trajectory, and it remains negligibly small for 
other bodies (as if we were considering 80 spacecrafts instead of one, but weakly coupled to 
the scalar field as in Brans-Dicke theory). 

In conclusion, although the above model (I5.10l) - fl5.17l) is not justified by any underlying 
symmetry principle, and should be considered more as a fit of Pioneer data than as a 
predictive theory, it proves that it is possible to account for this anomalous acceleration in a 
consistent field theory, while satisfying all other solar-system data (and binary-pulsar tests 
too, since the model reduces to Brans-Dicke theory in vacuum). Let us recall that there 
is no reason to trust Eq. (15.1 7p in a different regime than the one relevant to the Pioneer 
spacecrafts, and notably in cosmology. Note also that there exist many different ways to 
reproduce the same phenomenology, as exhibited by various expres sions of the d isformal 
coupling function B above. Our only addition with respect to Refs. 139l . 140l . 141 1 is that 
a varying post-Newtonian parameter 7 is indeed possible within a stable model admitting a 
well-posed Cauchy problem. 



VI. CONCLUSIONS 

One of the aims of this article was to clarify the mathematical consistency of various field 
theories proposed in the literature to reproduce the MOND phenomenology. We underlined 
that besides experimental constraints which need to be satisfied, including new ones that 
we discussed (notably binary-pulsar tests), there are also several basic requirements that a 
field theory needs to meet: stability (i.e., boundedness by below of its Hamiltonian) and 
well-posedness of its Cauchy problem (i.e., notably the fact that all field equations must be 



hyperbolic). In our opinion, the TeVeS model |30|, |31[ is presently the most promising, but 
it also presents several serious difficulties discussed in Sec. HVl and cannot be considered 
yet as a fully consistent theory. We examined another possible route to reproduce the 
MOND phenomenology, within the class of the already studied RAQUAL models [13], but 
with a different spirit. Its originality is that the theory is particularly simple in vacuum 
(pure general relativity or Brans-Dicke theory). However, the analysis of the field equations 
within matter exhibited a deadly inconsistency. Although this framework did not provide 
any serious alternative to TeVeS yet, we anyway believe further work along its line might do 
so. Its interest also lies in its relative simplicity. The model proposed in Sec. IV At assuming 
a nonminimal coupling of matter to the spacetime curvature, is for instance an excellent 
toy model of other MOND-like field theories, as its exhibits their generic difficulties: non- 
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predictiveness, fine-tuning, and above all instability. 

It should be underlined that other experimental problems also need to be addressed by 
any field theory of MOND, although we did not discuss them in the present article. First 
of all, weak lensing observations have already exhibited several evidences of dark m atter 
not located around clusters of baryon ic rn atter, notably the famous bullet cluster^^ jlSOl 
and the very recent dark matter ring |l51 |. The basic idea of the MOND dynamics seems 



therefore to be contradicted. However, several studies |l52 . 153 . 154 . 155l | have argued that 



such observations do not rule out MOND-like field theories, whose predictions can differ 
significantly from the original MOND proposal in non-spherical and dynamical situations. 
Another generic problem of models avoiding the dark matter hypothesis is to predict the 
right relative heights of the second and third acoustic peaks of the CMB spectrum. Indeed, 
in absence of coUisionless dark matter, their heights should decrease monotonically by Silk 
damping, whereas WMAP data js], 0] confirm that the second and thir d peaks have similar 



heights (consistent with the dark matter paradigm). Reference |156| analyzed the CMB 



spectrum within the TeVeS model, and found that it can fit experimental data provided 
their is a present density fli, ~ 0.17 of massive neutrinos. This number is not far from 
the needed dark matter density Qdm ~ 0.24 in the standard ACDM (cosmological constant 
plus cold dark matter) cosmological model, and illustrates that the existence of dark matter 
(as massive neutrinos) would anyway be needed in the TeVeS model, although the MOND 
dynamics was initially devised to avoid the dark matter hypothesis. However, it is already 
known that some amount of dark matter is anyway needed to account for the cluster rotation 
curves, within the MOND model, and massive neutrinos would be perfect candidates since 
they can cluster on such scales although they are too light to cluster on galaxy scales. 
Therefore, a MOND-like field theory remains a possibility, provided there is a large enough 
density of massive neutrinos — or any other dark matter candidate with similar properties. 
In the model [33|, closely related to TeVeS, one of the two scalar degrees of freedom is 
assumed to play the role of this dark matter. A third difficulty of MOND-like field theories 
is also to explain the numerical coincidence ~ cHq/G between the MOND acceleration 
constant and the Hubble constant (i.e., the present expansion rate of the Universe) (l7| . 
Some very promising ideas have been developed notably in [32,] , but it cannot be considered 
as a prediction. For instance, it may happen that Oq is actually related to the cosmological 
constant instead of the Hubble expansion rate, via a relation like ao ~ (-\/A)c^/8. In such a 
case, it would not vary with time, and one would "only" have to explain why it happens to 
take this precise value. On the other hand, a constant value of would probably be ruled 
out by the observed CMB spectrum, although this obviously depends on the precise features 
of the model. A fourth and related difficulty is to explain structure formation without cold 
dark matter, while remaining consistent with the tiny amplitude of CMB fiuctuations. As 
far as we are aware, there is at present no obvious solution to this crucial problem, and 
a larger value of at early cosmological times might be necessary. The modification of 
gravity, the possible renormalization of Newton's constant G at large distances, and the 
contribution of the various fields entering the model might also contribute to a solution. A 
fifth difficulty is more theoretical than experimental. The action-reaction principle implies 
that if a mass mi feels a force {y/Gm2ao)mi/r caused by a second mass m2, in the MOND 
regime, then the latter must feel the opposite force and thereby undergo an acceleration 



On the other hand, the cluster Abell 520 
within MOND than as dark matter [149 1 



148l | exhibits a dark core which might be more natural to explain 
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{■\/ GaQ/ m2)mi/ r . Since this expression diverges as m2 — > 0, any light enough particle should 
thus be infinitely accelerated by distant objects feeling its MOND potential. Actually, the 
self-energy of such particles must be properly renormalized, and it has been proven in [13] 
that the RAQUAL model does not suffer from any inconsistency due to the action-reaction 
principle. Indeed, the force felt by mass mi takes the MOND form i[_^/Gm2ao)"^i/''^ only 
when 1712 ^ ^1- The exact expression of the forces has been derived in 157l |. miOi = 771202 = 
(2v/Gao/3r)[(mi + m2)'^/^ — m^^^— mg'^^], so that the acceleration 02 = \/GmiaQ/r + 0{^/rri2) 
remains finite even when 1112 0. Although we expect that similar arguments would apply 
for any field theory deriving from an action principle, it remains to explicitly prove so for 
new models that one may consider, for instance for the class of nonminimally coupled fields 
studied in Sec. |V] above. 

All these difficulties, in addition to those discussed in the present paper, underline that 
the construction of a consistent field theory reproducing the MOND dynamics and the right 
light defiection is far from being obvious. Although the considered actions depend on several 
free functions allowing us to fit different kinds of data, there are so many theoretical and 
experimental constraints that no model passes all of them at present. The conclusion of 
our analysis seems thus to be in favor of the dark matter paradigm. On the other hand, 
the Tully-Fisher law, i.e., the existence of a universal acceleration scale ao, has not yet been 
derived in dark matter models. More generally, dark matter profiles seem to be tightly 
correlated to baryonic ones jl58l | (up to apparent counterexamples such as the bullet cluster 
and this lacks any explanation in CDM whereas this is an obvious prediction of 



MOND. Therefore, the dark and baryonic matter clustering still needs to be understood in 
more detail, and the analysis of modified gravity theories remains an interesting alternative 
to the standard ACDM paradigm. 

Identifying theoretical difficulties and finding ways to solve them often provides new 
possible interpretations of experimental data. For instance, we saw in Sec. Ill Fl that the 
RAQUAL kinetic term needed to be slightly modified for vanishingly small accelerations, in 
order for the field theory to remain consistent. The resulting model predicts that gravity 
becomes Newtonian again at very large distances (after a transition by the MOND regime), 
but with a much larger gravitational constant. Both the cosmological predictions and our 
interpretation of cosmological data would be fully changed within such a model. This may 
lead to a new understanding of both the "local" mass discrepancies (within galaxies and 
clusters) and of the cosmological dark matter. Another bonus of our study was also to 
provide an example of a consistent field theory reproducing the Pioneer anomaly without 
spoiling other predictions of general relativity. But although it passes theoretical and ex- 
perimental constraints, it does not pass "esthetical" ones (as defined in the Introduction), 
i.e., it is actually tuned to account for Pioneer data. It would remain to find an underlying 
symmetry principle which could predict this phenomenology. 
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